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tions in  simple  cubic  crystals  is  presented.   The  crystal  is  considered 
to  be  a  continuum  where  defects  are  simulated  by  arrays  of  point  forces 
acting  on  the  centers  of  atoms  in  the  immediate  neighborhood  of  the 
defect.   The  magnitude  of  these  forces  is  determined  by  the  condition 
that  they  have  the  same  displacement  field  as  the  corresponding  defect 
in  the  ordinary  continuum  model.   Infinitesimal  prismatic  and  shear 
loops  are  constructed  for  simple  cubic  crystals  and  used  to  construct 
screw  and  edge  dislocations.   The  arrangement  of  the  atoms  in  the 
vicinity  of  the  dislocation  line  is  obtained  and  compared  to  Peierls' 
model.   The  self -energy  of  these  dislocations  is  found  to  be  of  the 
correct  form  provided  the  force  cbnstants  are  correctly  determined. 

Atomic  arrangements  around  kinks  in  screw  and  edge  dislocations 
have  been  computed  and  are  presented.   The  model  developed  promises  to 
be  of  great  value  in  studying  atomic  displacements  in  the  vicinity  of 
the  dislocation. 
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CHAPTER  1 


INTRODUCTION 


A  thorough  comprehension  of  the  nature  of  defects  in  crystalline 
materials  and  especially  in  metals  is  fundamental  in  order  to  explain 
many  of  the  properties  and  the  behavior  of  these  solids.   In  particular, 
vacancies,  interstitials  and  dislocations  cannot  be  ignored  when  diffu- 
sion, mechanical  behavior,  electrical,  optical  and  magnetic  properties 
are  studied. 

The  usual  theory  of  lattice  defects  assumes  a  "local"  continuum 
model.   The  matter  concentrated  in  the  atoms  is  supposed  to  be  uniformly 
distributed  over  the  whole  space  occupied  by  the  crystal.   The  local 
atomic  arrangement  is  ignored  and  the  defect  is  replaced  by  a  singular 
line,  point  or  surface  in  a  continuum  body  [1-3J. 

This  model  has  proved  to  be  extremely  valuable  for  studying 
properties  which  are  not  sensitive  to  the  atomic  configuration  in  the 
vicinity  of  the  defect,  but  it  is  limited  by  the  discrete  nature  of  the 
atomic  array.   Consequently,  it  is  always  necessary  that  expressions 
for  the  displacement  field  of  the  defect  be  terminated  at  some  distance 
from  it.   Furthermore,  since  the  continuum  approximation  ignores  the 
local  atomic  arrangement  around  defects,  it  disregards  the  short  range 
anisotropy  of  the  displacement  field. 

To  remedy  these  shortcomings,  atomistic  computations  have  been 
attempted.   Atomic  positions  and  interactions  are  considered  explicitly 


in  the  core  region  of  the  defect,  with  some  laws  defining  the  pairwise 
atomic  potentials  (see  [4-13]).   Further  from  the  center  of  the  defect, 
continuum  theory  is  assumed  to  hold,  so  that  the  only  atoms  which  need 
to  be  considered  are  those  whose  positions  are  necessary  for  calcula- 
tions of  energies  for  the  core  region.   Such  a  method  involves  first 
the  construction  of  a  suitable  interatomic  potential,  and,  secondly, 
sums  over  a  large  number  of  lattice  points  which  have  to  be  carried  out 
numerically.   It  is  undoubtedly  the  best  existing  method  of  determining 
the  local  atomic  arrangement  around  defects,  but  it  is  costly  and  very 
sensitive  to  the  chosen  size  of  the  core  region  [14].   Furthermore,  it 
involves  convergence  problems,  and  the  manner  in  which  boundary  condi- 
tions are  imposed  is  very  delicate. 

It  is  therefore  worth  exploring  methods  refining  the  ordinary 
continuum  model  by  introducing  the  atomic  arrangement  of  the  crystal, 
but  with  a  minimum  increase  in  computational  effort.   In  such  a  model, 
the  atoms  will  be  considered  to  be  embedded  in  a  continuum  and  the 
defect  formed  by  the  placement  of  suitable  point  forces  at  positions 
corresponding  to  atomic  sites  close  to  the  defect  [15].   The  result- 
ing displacement  field  is  the  sum  of  the  displacement  fields  of  all 
the  point  forces  and  is  taken  as  the  displacement  field  around  the 
defect.   Examples  of  such  constructions  by  superposition  of  infini- 
tesimal loops  have  been  given  by  Koehler  [16],  Groves  and  Bacon  [17] 
and  Kroupa  [18]  for  local  continuum  models. 

In  this  present  study,  we  shall  concern  ourselves  in  examining 
straight  screw  and  edge  dislocations  in  simple  cubic  crystals.  First, 
a  brief  description  of  point  forces  and  infinitesimal  primitive  loops 


will  permit  us  to  analyze  the  displacement  field  and  self-energy  of  dis- 
locations.  The  screw  dislocation  will  be  constructed  from  an  array  of 
primitive  shear  loops,  and  the  edge  dislocation  from  an  array  of  prim- 
itive shear  loops  and  prismatic  loops.   Comparison  will  be  made  between 
both  models  in  the  case  of  the  edge  dislocation.   Furthermore,  calcula- 
tions of  atomic  displacements  around  kinks  will  be  attempted  for  both 
dislocations. 


CHAPTER  2 
BASIC  CONCEPTS 

Point  Force 

A  point  force  F  is  a  highly  localized  body  force  distribution 
applied  to  a  material  point  in  a  continuum. 

F.(r')  =  f.  f  6(R)  dV-  (1) 

J     J  v  r 

where  R  =  r  -  r'  and  6 (R)  is  the  Dirac  delta  function.   The  displace- 
ment field  at  r,  u  (r) ,  due  to  such  a  point  force  at  r'  ,  can  be 
i 

obtained  from  the  equilibrium  equations  of  elasticity  and  Hooke1 s  law 
[2]  in  the  form 

u.(r)  =  F.(r',)G.  .  (R)  (2) 

J       i     ij 

where  for  an  infinite  isotropic  body 


G..(R)=       * 


ij  16rrp, 


^    |r|  +  ^  |r|3_ 


(3) 


G   (R)  is  the  Green1  s  tensor  response  function  for  a  point  force, 
ij 

It  is  the  component  parallel  to  x .  of  the  displacement  field  at  r'  due 

to  a  unit  point  force  parallel  to  x.  at  r.  It  can  be  shown  that  the 

Green' s  tensor  is  symmetric. 

As  we  see  from  the  expression  for  G  . (R) ,  this  function  is  not 

U 

defined  for  R  =  0  and  we  are  unable  to  determine  the  displacement  of 


the  point  of  application  of  the  point  force  from  Equation  (2).   To 
remove  this  mathematical  divergence,  we  shall  associate  a  finite  dis- 
placement u(r')  with  the  point  of  application  of  the  point  force. 
This  value  u(r')  can  be  considered  as  being  the  average  resultant  dis- 
placement of  points  on  a  surface  surrounding  the  point  of  application 
of  the  point  force.   This  surface  is  determined  such  that  u(r' )  is  the 
mean  value  of  the  vector  displacements  of  two  points  symmetrical  with 
respect  to  r' .   This  vector  is  acting  in  the  same  direction  as  the 
point  force. 

u(r')  =  fpCr'  +  r )  +  u(r'  -  rY]  .  (4) 

Using  polar  coordinates  as  shown  in  Figure  1,  the  absolute  value  of  the 
displacement  u(r' )  takes  the  following  form. 

2„-l 


|u(?')|  = 


16TTH 


3-4V  1     1   cos 

—  + 


(5) 


l-v   r    l-v   r 

So,  for  a  definite  value  of  ju(r')]  we  can  define  a  surface  of  revolu- 
tion about  the  direction  of  F,  on  which  all  the  points  have  displace- 
ment components  |u(r')|  along  F.   This  surface  surrounds  a  volume  which 
can  be  considered  as  a  core  surrounding  the  point  of  application  of  the 
point  force.   The  core  can  be  interpreted  as  the  volume  where  Equation 
(2)  for  the  displacement  field  is  r.o  longer  valid.   The  size  of  the 
core  depends  directly  on  the  value  of  Ju(r')J  assigned. 

It  must  be  pointed  out  that  the  average  value  of  the  radius 
vector,   rj ,  of  the  core  is  equal  to  the  radius  r  of  the  sphere  on 
which  the  average  displacement  of  its  points  is  equal  to  Ju(r')|, 
that  is 


°   24TT^|u(r')| 

On  the  other  hand,  |u(r')|  can  be  related  to  a  force  constant.   By  anal- 
ogy with  a  discrete  lattice  model,  such  a  point  force  applied  on  an 
atom  causes  an  equal  and  opposite  resisting  force  proportional  to  the 
displacement  of  the  atom  given  by 

|f|  =  c|u(r*')|  .  (7) 

C  is  known  as  a  force  constant,  and  is  the  force  acting  on  an  atom 
required  to  produce  a  unit  displacement.   In  other  terms,  its  inverse 
is  the  displacement  of  the  atom  caused  by  a  unit  force  acting  on  it. 

This  force  constant  is  the  parameter  we  shall  use  in  the  follow- 
ing problems  encountered.   It  will  be  determined  for  each  special  case 
by  requiring  that  our  mathematical  model  obeys  certain  physical  imper- 
atives.  It  will  be  straightforward  to  deduce  |u(r')|  and  the  size  of 
the  core  from  the  value  of  C. 

The  self-energy  of  a  point  force  is  defined  as  being  the  work 
done  by  this  force  against  interatomic  reaction  forces  when  it  is 
introduced  into  the  continuum.   So,  using  Equation  (7) 

Es  =  |f|u(?')|  =|L..  (8) 

The  interaction  energy  between  two  point  forces  F    (r1 )  and 

FC2)(r')  is 

m 

ET  =  -  F,(1)(r)  u(2)(r')  =  -  F(1)  (?)  F(2)(r*')  G   (?-?'),        (9) 

I      k       k  k        m        km 


where  the  sign  is  determined  following  Cottrell's  convention  for  dis- 
location interaction  energies,  i.e. ,  it  is  the  work  done  by  external 
forces  when  the  second  force  is  applied  in  the  presence  of  the  first, 
or  vice  versa.  The   total  elastic  potential  energy  of  the  system  is 
the  sum  of  the  self-energies  of  the  two  point  forces  less  the  pairwise 
interaction  between  them. 

Double  Force 

A  double  force  is  constructed  from  two  eqiial  and  opposite  point 
forces  ±  F  applied  at  points  separated  by  a  vector  h.   If  the  forces 
are  collinear,  we  have  a  double  force  without  moment,  otherwise  with 
moment.   The  strength  of  a  double  force  is  defined  as 

P*k=   lira     (hiV-  (10) 

|h|-o 

|f|  — 

The  displacement  field  is  obtained  by  superposition.   If  the 
separation  distance  between  the  forces  is  very  small,  we  can  expand 
the  displacement  field  of  each  in  a  Taylor  series  about  the  midpoint 
of  h.   Keeping  only  the  first  order  terms,  we  are  led  to  the  displace- 
ment field  mentioned  by  Kroner  [2] 

u.(7)  =P£k(?<)  Gk.y ?-?•).  (id 

As  for  a  single  point  force,  the  displacement  of  the  point  of 
application  is  undefined,  but  this  divergence  can  be  removed  in  the 
same  way  as  before  by  introducing  the  concept  of  a  core  surrounding 
a  double  force. 


Primitive  Dislocation  Loops 

Following  Kroner's  definition  [2],  an  infinitesimal  dislocation 
loop  in  a  continuum  is  the  boundary  of  a  microscopic  surface  which 
separates  regions  in  the  continuum  which  have  suffered  a  relative  dis- 
placement b.   The  Burgers  vector  of  the  loop  is  defined  as  the  line 
integral  of  the  elastic  displacement  u  around  a  circuit  containing  the 
dislocation.   The  displacement  field,  at  a  point  r,  of  such  a  loop  of 
surface  dS  with  normal  vector  n  and  centered  at  r'  is  found  to  be 

n  (r)  «  b.n.  c...  .  C.  .   dS  (12) 

m       1  j   ijk£  km,x 

where  c.  .,  .  are  the  elastic  constants, 
ljki 

The  similarity  between  this  expression  and  the  displacement 
field  of  a  double  force  (Equation  (11))  leads  us  to  consider  the  infin- 
itesimal dislocation  loop  as  a  nucleus  of  strain  with  the  fundamental 
double  force  tensor 

P.  #  =  c  ...  b.n.  dS  (13) 

kX    ljkX  i  j 

or,  for  an  isotropic  continuum, 

\i  -  CXVji  +  6i/j*)  +  xVk  J  Vj  ds  (14) 

Up  to  this  point,  we  have  completely  ignored  the  local  atomic 
arrangement  around  the  loop.   In  real  crystals  the  interatomic  reaction 
forces,  developed  when  the  atoms  are  displaced  to  form  the  defect,  are 
the  physical  origin  of  the  double  force  tensor  characteristic  of  the 
dislocation  loop.   So  it  seems  logical  to  construct  such  a  loop  by 
applying  point  forces  in  the  continuum,  but  at  points  corresponding 
to  atomic  positions  located  immediately  around  the  defect.   The 


displacement  field  of  the  loop  is  then  the  superposition  of  the  dis- 
placement field  of  each  point  force.   Each  primitive  loop  has  the 
character  of  a  "unit  cell"  for  the  defect.   These  "unit  cells"  can  be 
assembled  to  form  a  more  complicated  defect  like  dilatation  centers  or 
dislocations.   So  logically  we  can  characterize  the  surface  of  the  loop 
dS  such  that  the  produce  |b  •  dS|  equals  one  atomic  volume  in  the  crystal 
structure  considered.   This  procedure  will  be  analyzed  more  specifically 
for  simple  cubic  crystals. 

(a)  Primitive  prismatic  loop  in  simple  cubic  crystals 

The  arrangement  of  the  first  neighbors  of  a  vacancy  loop  in 
simple  cubic  crystal  is  shown  in  Figure  2.   A  primitive  prismatic  loop 
is  constructed  in  the  following  steps.   First  a  vacancy  is  created  by 
removing  an  atom  from  the  lattice.   This  vacancy  is  simulated  by  apply- 
ing on  its  first  neighbors  forces  of  magnitude  F  directed  towards  the 
vacancy  center.   In  the  second  step,  two  extra  forces,  ±  G,  are  applied 
in  a  direction  normal  to  the  {00l}  plane,  on  the  atoms  in  the  (OOl) 
direction,  towards  the  center  vacancy  in  order  to  collapse  the  configur- 
ation onto  the  {001}  plane.   In  this  manner,  we  have  set  up  three 
double  forces,  all  without  moment,  leading  to  the  diagonal  dipole 
tensor 

Pll  =  P22  =  2Fa' 
and  (15) 

P33  =  2a(F+  G), 

where  a  is  the  lattice  parameter  of  the  simple  cubic  crystal. 

Though  each  pair  of  forces  is  clearly  separated  by  a  distance 
2a,  at  distances  from  the  loop  large  with  respect  to  the  interatomic 
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Figure  1.   Coordinate  System  to  Evaluate  the  Core  Region 
Around  a  Point  Force 
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Figure  2.   Prismatic  Loop  in  Simple  Cubic  Crystal 


11 


distance  they  appear  as  three  double  forces  which  can  be  identified 
with  a  dislocation  loop  as  described  above.   The  Burgers  vector  of 
this  loop  must  represent  the  collapse  of  the  atoms  in  the  (OOl) 
direction,  whose  relative  displacement  must  be  a  in  order  to  create 
a  new  regular  arrangement  of  the  atomic  planes.   As  stated  previously, 
the  surface  dS  is  chosen  such  that  |b  •  dS|  is  equal  to  a  here.   So, 
following  Equation  (14) 

Pll  =  P22  =  X*3 
and  (16) 

P33  =  a  +  2^}  a  ' 

The  forces  applied  on  the  atoms  can  now  be  obtained  by  comparing 
Equations  (15)  and  (16).   The  displacement  field  and  the  self-energy 
of  the  loop  can  easily  be  deduced. 

(b)   Shear  loop  in  simple  cubic  crystal 

A  primitive  shear  loop  in  a  simple  cubic  lattice  (Figure  3) 
is  constructed  as  follows.   The  forces  F  applied  on  the  atoms  impose 
the  direction  of  the  shear.   Since  the  loop  must  be  kept  in  equilib- 
rium with  respect  to  its  center,  additional  forces  G  have  to  be  applied, 
forming  a  couple  whose  moment  about  the  center  of  the  loop  counter- 
balances that  of  the  shear  forces  F.   The  Burgers  vector  of  the  loop 
is  the  smallest  shift  allowed  by  the  atomic  arrangement.   By  the  same 
method  as  for  the  prismatic  loop,  it  is  found  that  the  only  nonvanish- 
ing  components  of  the  dipole  tensor  are 

P12  =  P21  =  ^  =  2Fa  =  2Ga'  (17) 
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so  the  magnitude  of  the  forces  has  the  value 

2 
F  =  G  =  p,  ~   .  (18) 

The  shear  loop  can  be  represented  with  respect  to  its  principal  axes 

x'  and  x'  (Figure  4)  leading  to  the  dipole  tensor 
1      ^ 

Pil  =  2aF  =  "  P22  •  (19) 

which  represents  two  double  forces  without  moment,  perpendicular  to 
each  other  and  acting  in  opposite  senses, 
(c)   Conclusion 

The  primitive  dislocation  loops  as  described  above  are  the 
basic  elements  for  our  process  of  simulating  larger  defects,  especi- 
ally dislocations.   We  shall  see  that  a  suitable  array  of  shear  loops 
can  describe  either  an  edge  or  a  screw  dislocation,  but  that  an  array 
of  prismatic  loops  can  only  simulate  an  edge  dislocation. 
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Figure  3.   Shear  Loop  in  Simple  Cubic  Crystal 
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Figure  4.   Shear  Loop  with  Principal  Axes 


CHAPTER  3 


RECTANGULAR  DISLOCATION  LOOP  IN 
SIMPLE  CUBIC  CRYSTAL 


A  rectangular  dislocation  loop  having  a  Burgers  vector  a(lOO) 
can  be  simulated  by  a  rectangular  array  of  primitive  shear  loops, 
stacked  as  shown  in  Figure  5.   The  dimensions  of  the  array  are  con- 
sidered to  be  very  large  compared  to  the  atomic  distance.   The  axes 
of  reference  are  shown  in  Figure  5  with  their  origin  at  the  center  of 
the  loop.   In  this  chapter,  we  are  only  interested  in  obtaining  prop- 
erties of  the  rectangular  dislocation  loop  related  to  our  main  inter- 
est, the  displacement  field  and  self -energy  of  the  pure  screw  and  edge 
dislocations. 

Displacement  Field 

The  displacement  field  at  any  points  of  this  array  is  simply  the 
sum  of  the  displacement  field  of  each  point  force 

u  (?)  =   S   F.(r')  G.  (r-r* ),  (20) 

m      .  -\   i      1m 

where  G   (R)  is  defined  in  Equation  (3).   Developing  this  sum  leads  to 
im 

the  general  expression 


2   R/2   I  /  L/2      L/2-1 
E   /  [   S   +     E 

p=-R/2|  lq=-L/2    q=-(L/2-l) 


u  (r)  =  u.  -g-  Z       (Is    +      z 


G2m(xrpa'X2-qa' 
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G3m(Xrpa'    V¥'    X3-|}    +    G3m(Vpa'    X2  +  ¥'    X3  +  f  > 


'    lG3m(Xrpa>    X2"¥'    X3-f}    +    G3m(VPa'    X2"¥'    X3+l} 


(21) 


These  equations  for  each  component  are  valid  everywhere  in  the  con- 
tinuum.  They  can  be  simplified  for  each  particular  region  of  the  loop. 

The  regions  of  the  continuum  where  a  pure  screw  dislocation 
is  simulated  correspond  to 


xi  =  "  T  +  ei 


X2    2 


and 


Ra 
Xl  =  T  ~   €i 


X2=  62 


where  e   and  e   are  small    compared    to  Ra  and  La,  respectively. 
Similarly,  regions  where  the  loops  have  a  pure  edge  character  corre- 
spond to 


and 


Xl=  ei 


La 
X2  =  T  "  62 


where  e   and  e   are  small   compared    to  Ra  and  La,  respectively. 
Each  particular  case  will  be  considered  in  the  following  chapters. 


Elastic  Potential  Energy 

The  work  done  by  the  forces  comprising  the  array,  that  is,  the 

energy  of  the  system,  is  the  sum  of  the  self-energy  of  each  point  force, 

minus  the  pairwise  interaction  energies,  as  defined  by  Equations  (8) 

and  (9).   We  call  W    the  energy  of  a  row  of  forces,  that  is,  two  lines 
Row 
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of  forces  parallel  to  the  x   axis  for  a  given  x   coordinate,  and  we  call 

W    the  interaction  energy  between  two  rows  of  forces,  as  defined  above. 
Int 

Following  these  appellations,  the  total  energy  of  the  system  has  the  form 


R 

W  =  (R+l)Wn   +   E   (R+1-P)WT 

Row     „         Int 
P=l 


(22) 


W    and  W    have  the  following  expressions,  respectively, 


3  5 
(J,  a 


Row 


2L-1 


L  C2     C3J 


2  4 


+  1Sr-CG33(0'0'a)  "  (2L-1}  G22(°'°'aD 


+  Ji~-  (g   (O.La.O)  -  G22(0, La, a)  -  G33(0,La,0)  -  G33(0,La,a)] 
2  4  L      r  -i 

+  2(13    E  (L-q)  [_G22(0,qa,0)  -  G22(0,qa,a)J 


q=l 

2  4  2  4 

-  n  a  G23(0,La,a)  -  2^  a   E  G23(0,qa,a) 

q=l 


(23) 


Int 


^2a4(2L-l)  (G22(pa,0,0)  -  G22(pa,0,a)]  +  p,2a4  Q}33(pa,0  ,0)+G33(pa,0,a)] 
+  u-  a   Q}22(pa,La,0)  -  G22(pa,La,a)  -  G33(pa,La,0)  -  G33(pa,La,a)J 


2  4  2  4 

2p,  a  G0„(pa,La,a)  -  4p,  a   E  G   (pa,qa,a) 
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L 


Q=l 
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2  4         r~  ~\ 

+  4|j,  a   E  (L-q)  [G      ipa,qa,0)    -   G22(pa,qa,a)J  . 

q=l 
Following  Equation  (22)  ,  W  becomes 

W=  n,2a4(R+l)L^  -  G22(0,0,a)f2   E  [g^O, qa,0)    -  0^(0, qa, a)] 


(24) 


q=i 


R   L 


+  2   E  [G22(pa,0,0)  -G22(pa,0,a2)+4   E   E  (jJ^Cpa.qa.O) 
p=l  "        '     p=l  q=l 

"  G22(pa'qa'a[l}  +  ^4^D  V?(^~C^)+  l&22(0'°'a)  +G33(°'°'a>] 

+  §  (G22(0,La,0)  -G22(0,La,a)  -G33(0,La,0)  -  G33(0,La,  a)] 


2   E  q  [p22(0,qa,0)  -  G29(0,qa,af| 
Q=l 


18 
R 


-   £  (G22(pa,0,0)  -G22(pa,0,a)  +  G   (pa.0,0)  +G  3(pa,0fa)) 
P=l 

R 
+   £  QG22(pa,La,0)  -  G22(pa,La,a)  -  G^Cpa  ,La,0)  -G   (pa.La.af) 


P=l 

R   L 


-  4   E   E  q  QG22(pa,qa,0)  -  G22(pa,  qa,  a)J  -G23(0,La,a) 

p=l  q=l 

L  R+l  R+l   L 

-2   E  G23(0,qa,a)  -2   £  G23(pa  ,La,  a)  -  4   £   E  G23(pa,qa,a) 
q=l   "  *  p=l  p=l  q=l 

+  ^a4L<-2  Ep  [G22(pa,0,0)  -G22(pa,0,a)] 

L   p=l 

R   L  "\ 

-4   E   E  p  [p22(pa,qa,0)  -  G^pa.qa, af)  > 
P=l  q=l  "'  J 

+  n  a  <   E  p[G22(pa,0,0)  -G22(pa,0,a)  -G33(pa,0,0)  -G33(pa,0,a)] 

-  E  pQ}22(pa,La,0)  -G22(pa,La,a)  -G33(pa,La,0)  -G33(pa,La,a)J 


P=l 

R   L 
+  4   E   E  pq  [_G22(pa,qa,0)  -  G22(pa,qa,a)J 
p=l  q=l 

R+l  L   R+l 

+  2   E  p  G23(pa,La,a)  +  4   E   E  p  G23(pa,qa,a) )   .       (25) 
p=l     "  q=l  p=l 

The  only  mathematical  difficulty  in  the  computation  of  such  a  formidable 

expression  lies  in  computing  the  single  and  double  sums.   The  first  terms 

are  computed,  up  to  a  chosen  integer  N  (usually  N=20),  and  the  rest  of  the 

terms  are  approximated  by  an  integral.   The  following  approximation  follows: 

L  L 

E  f  (q)  «  f(q)  +|  f  (N)  +|  f(L)  +J  f  (x)  dx,  (26) 

q=l  "        ""  N 

and  for  a  double  sum: 

R   L  N-l  N-l  N-l 

E   E  f(p,q)  m  E   E  f(p,q)+|  E  Qf(p,N)  +  f(p,L)) 
p=l  q-=l  p=l  q=l  p=l 

N-l  _  N-l   L 


+  |  E   ff(N,q)  +  e(R,qf|  +   E  J  f  (p,y)  dy 
q=l  ^~  p=l  N 
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+   E  J   f (x,q)  dx  +  i  f  f (N,y)  dy  +  i  J   f (x,N)  dx 

q=l  N  N  N 

-   L  R  .. 

+  ±  f  f  (R,y)  dy  +  i  j  f  (x,L)  dx  +  -  [_f  (N,L)  +  f  (R,N) 
N  N 

L  R 

+  f (N,N)  +  f (R,L)j  +  J  J  i (x,y)  dx  dy  . 


(27) 


These  approximations  give  the  correct  form  of  the  divergence 

4 

in  R  and  L  for  divergent  sums,  and  give  an  accuracy  of  1  part  in  10 

when  the  sum  converges,  which  is  sufficient  for  the  model  employed. 


The  final  expression  for  W/ua  becomes 
(R+1)L   .   ua   2.8545  -  2.3676V      R+l 


pa 


L 

+    4- 


4tt 


„       ua      2.  8545    -   2.3676V 

4tt  J— ■  + 

C2  1    -    V 


2tt(1-v) 


2n 


2RL 


2n 


C~ 


+  «/R  +l 


R+l 

4tt 


R 


+  ^ 


2      2 

+L 


.   1         1  \         1.2979    -.6458V 


1.6420  +  .3580V 


1    -    V 


2(2-y)   Jr  +L 
1  -  V  4tt 


/X    2 


'        fl         1  A       3.3903  -  .6872V 
2ttu^---)  - 


4tt(1-v) 


(28) 


W  is  not  the  energy  of  the  physical  dislocation  loop.   It  con- 
tains an  extra  strain  energy  in  the  region  bounded  by  the  plane  on  which 
forces  are  applied  and  the  boundaries  of  the  array,  which  does  not 
account  for  the  relaxation  of  the  atoms.   In  the  way  the  forces  have  been 
applied,  a  relative  displacement  greater  than  a/2  has  been  created  across 
most  of  the  slip  plane  for  the  atoms  reaching  their  final  configuration. 
It  is  from  this  final  configuration  that  the  relative  displacement  proce- 
dure must  be  measured  in  order  to  calculate  the  actual  strain  energy 
stored  between  the  planes  of  forces.   Such  a  correction  energy  will  be 
computed  in  each  case,  for  the  pure  screw  and  pure  edge  dislocation. 

We  shall  not  examine  all  the  properties  of  the  dislocation  loop 
here,  since  our  purpose  is  to  treat  this  loop  as  an  intermediate  step 
in  simulating  screw  and  edge  dislocations. 


CHAPTER  4 
SCREW  DISLOCATION  IN  SIMPLE  CUBIC  CRYSTAL 

From  the  results  found  in  the  previous  chapter,  the  displacement 
field,  the  relative  displacement  across  the  slip  plane,  and  the  self- 
energy  of  the  screw  dislocation  in  simple  cubic  crystal  will  be  obtained. 

Displacement  Field 

As  seen  previously,  the  regions  of  the  dislocation  loop  having 
the  characteristics  of  a  pure  screw  dislocation  correspond  to 


r 


Ra 

T+   ei 


Ra 

xi  =  T  ~   ei 


and   < 


2    2 


V 


where  6   and  e   are  small  compared  to  Ra  and  La,  respectively. 

Since  both  regions  represent  two  identical  dislocations,  but  of 
opposite  signs,  we  shall  only  consider  the  first  one,  and  translate  the 
x  axis  by  an  amount  of  -Ra/2  so  that  it  becomes  the  boundary  of  the 
array.   We  shall  keep  the  same  symbols  x   and  x  for  the  new  variables. 
(See  Figure  6. ) 

The  only  change  in  the  expression  for  the  displacement  field  as 
written  in  Equation  (21)  is  that  p  is  now  summed  from  zero  to  R.  Since 
an  analytical  expression  is  desired  for  u  ,  an  approximation  different 
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21 
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from  that  given  by  Equation  (26)  will  be  employed  to  compute  the  dis- 
crete sums.   Euler's  formula  [191    is  most  suitable  for  this  case: 

b         b 

E  f  (P)  =  J   f  (x)  dx  +  -  (_f  (a)  +  f  (bf]  +  —  (V  (b)  -  f '  (a)J 

a        a 

.^[f...(b).f...(a)]+s_|_Qf«)(b).^u)3 


+  ...  .  (29) 

The  accuracy  of  the  approximation  depends  on  the  number  of  terms  used. 
The  advantage  of  this  method  is  that  at  most  of  the  points  where  the 

displacement  field  is  computed,  the  three  first  terms  are  sufficient 

-3 

for  the  accuracy  required,  i.e.,  a  relative  error  of  10   is  accepted. 

The  first  summations  which  have  to  be  computed  are  the  two 
summations  on  q.   These  sums  can  be  computed  exactly  because  we  always 
consider  x   small  with  respect  to  La,  or  in  other  words,  our  range  of 
interest  is  far  from  both  ends  of  the  dislocation  line.   So  we  shall 

have 

L/2  L/2-1  +  °° 

2     f(q)  =      S       f  (q)  =  J        f  (x)  dx  (30) 

q=-L/2         q=-(L/2+l)         -« 


where  f  represents  the  whole  expression  to  be  summed.   The  following 
components  of  the  displacement  field  have  been  found: 

r 


WW  =  ° 


<  WW  =  4W  \^ 


oW  w  =  ° 


(x  -pa)   +  (x3+2-) 

2        a~~2 
p-0    (xi-pa)   +  (x3--) 


R 
a   „ 


(31) 
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We  can  already  notice  that  some  of  the  characteristics  of  the  displace- 
ment field  around  a  screw  dislocation  are  displayed,  that  is,  the 
components  u,  and  u„  are  equal  to  zero,  and  u   is  independent  of  x  . 

Euler' s  formula  can  be  applied  a  second  time  to  obtain  a  final 
expression  for  u  ,  the  accuracy  depending  on  the  number  of  terms  retained. 
We  shall  label  the  different  terms  composing  u   in  the  following  way, 


u2(xi'W  =  ^{g(Vpa)  dp  (32) 

vwv  =  £?  [g(xi}  +  g(xrRaD  (33) 

U2(X1' W  =  48^  0  (VRa)  -  g'  (X1}]  (34) 

U23(X1'X2'X3)  -  "  28^  &(3)<VR3)  "  ^^V]  '    (35) 


etc. ,  where 


2        a  2 
(x1-pa)   +  (Xg  +  g) 

g^-pa)  =  Sm g ^~2    '  ^36^ 

(xrpa)   +  (x3--) 

The  computation  of  u   raises  a  mathematical  problem,  since  we 
are  integrating  over  a  region  where  the  integrand  contains  singular 
points  for  certain  values  of  x  and  x   (x  =  na,  x  =  ±  a/2).   For  these 
points,  the  integral  can  be  broken  into  two  parts: 


Ra 
J  g(na-pa)  dp  -   lim 
0  e  -0 

For  this  specific  case 


"  na-e  Ra 

J  g(na-pa)  dp  +  j     g(na-pa)  dp 
0  na+e 


(37) 


2     2 
/     ^   n     (na-pa)   +  a  foR-j 

g(na-pa)  =  m  g —  •  ( 

(na  -  pa) 
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Since 


J  0n   (na-pa)   dp  =  -2  f(n-p)  Ski   (na-pa)  -  (n-p)~J  ,  (39) 

the  integral  is  equal  to  zero  for  p  =  n,  and  so 


Ra 


J   g(na-pa)  dp  =  G(na-Ra)  -  G(na). 
0 


(40) 


The   integration  makes    all    the   singular  points   vanish,    except 
for   those   at   x     =0   and  x     =   ±  a/2.      So   the  final    analytic   expressions 


for  u     become 


U2    <VX3>    =    4^ 


2n^ 


.      -1  1  .-11 

+   tan       +   tan       


3      2 


(41) 


+    2x, 


ft      -1       Xl 
tan        

L  X3  +  J 


-   tan 


2         ,  a  2 

Xl  +    <x3+2>     "\ 

ci fe  "2 — ; — T2  >. 

X,     +     (X„   -75-)       J 


3     2' 


2         ,  a,2 

1    ,  ,  *    0     Xl+    U3  +  2} 

u2  (xi'V  =  i^ ^-§ — ; — rs  ■ 

Xl+    (X3-2} 


(42) 


U2    (X1'X3}    =    "   24^ 


r 


2  as2  2  ax2 

Xl   +    (X3-2}  Xl   +    (X3  +  2)    J 


(43) 


and 


VW    =   "  730TT 


4  KC«*b-5»  -4D    *i[3"vf>  -O 


(44) 


Then,    since  u    (x     x   )    is    the   sum  of   u    (x    ,x   )    for   every    i, 

£t        X        o  *^        X        o 

U2(X1'X3)    =   U2(X1'X3}    +    "2^1'V    +    U2(X1'X3} 


+   u2(xllX3)    + 


(45) 
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Equation  (41)  takes  three  different  forms,  following  the  region  where 
x   is  computed:   rr  corresponds  to  x  >  a/2,  2tt(x  /a)  to  -a/2  <  x  <  a/2, 
-TT  to  x  <  -a/2. 

An  asyn.ptotic  expression  for  u   can  be  obtained  when  x   and  x 
are  considered  large  with  respect  to  the  atomic  distance,  but  still  far 
from  the  ends  of  the  dislocation  line. 


WV  =  2^ 


,  TT     ,-11 

±  —  +  tan   — 
-  2         X3J 


(46) 


This  is  the  well-known  expression  obtained  from  the  Volterra  solution 
for  a  screw  dislocation  [20] . 

The  relative  displacement  across  the  slip  plane  is  defined  as 


AW  =Vxi'2>  -  w-^  =  2Vxi'2> 


(47) 


It   can  be  directly   deduced  from  Equation    (31)    or   Equation   (42) 
to    (45): 


Au0(x  )   =      ±_       S    Sm 
21  2n     n_ 


2  2 

(x  -pa)      +    a 


p=0  (x1~pa) 


(48) 


0  12  3 

Au    (x  )    =   Au   (x  )  +.Au2(x  )  +  Au    (x   )  +  Au    (x  )   + 


(49) 


with 


AU2(X1}    =    2^ 


2  2 

,    x  x  x      +    a 

-11  1    n  1 

2    tan +   rr  +   —  m  = 

a  a  2 

x 
1        -J 


(50) 


2  2 

1                    a         Xl  +    3 
AU2(X1}    =    4^ 2— 


(51) 
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AVV  =  m 


2     2 
x  +  a 


1  J 


(52) 


and 


AU2(X1}  =  36Cm 


,  2 
xi(x1 


3a2)  n 


3     /  2     2\3 
Cl     (Xl  +  3  ) 


(53) 


The  more  and  more  precise  expressions  for  Au  (x  )  are  plotted 
in  Figures  7  and  8,  and  listed  partially  in  Table  1.   A  remarkable 
precision  is  obtained  for  the  regions  where  |x  |  ^  a   after  evaluating 
only  a  few  terms.   On  another  hand,  the  only  term  which  is  not  singular 
at  x  =  0  is  Au  (x  ).   It  seems  to  deviate  significantly  from  the  correct 
curve  for  Au  (x  ).   However,  the  order  of  magnitude  of  the  real  relative 
displacement  at  x   can  be  obtained  approximately  by  interpolation.   So 
the  relative  displacement  of  atoms  above  and  below  the  slip  plane  is 
known  everywhere  except  at  x  =  0. 


TABLE  1.   Relative  Displacement  of  Atoms  Across  the  Slip  Plane  for 
a  Screw  Dislocation 


x  /a 


Au°/a 


+  Au  /a 


+  Au  /a 


3 

+  Au2/a 


0000 


.0000 


5.0 

.0316 

.0347 

4.0 

.0394 

.0442 

3.0 

.0521 

.0605 

2.0 

.0766 

.0943 

1.0 

.1397 

.1948 

.5 

.2243 

.3524 

.5 

.7757 

.9037 

1.0 

.8603 

.9155 

2.0 

.9234 

.9412 

3.0 

.9479 

.9563 

4.0 

.9606 

.9654 

5.0 

.9684 

.9715 

.0000 

.0349 
.0446 
.0614 
.0970 
.2081 
.3949 

.8613 
.9022 
.9385 
.9554 
.9651 
.9713 


.0000 

.0349 
.0446 
.0613 
.0968 
.2074 
.3890 

.8696 
.9033 
.9387 
.9554 
.9651 
.9713 


1.0000 


1.0000 


1.0000 


1.0000 
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Figure  7.   Relative  Displacement  Near  the  Core  of 
a  Screw  Dislocation 
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Another  way  of  computing  the  relative  displacement  at  atomic 
points  is  to  go  back  to  the  definition  of  the  displacement  field  by 
Green' s  functions. 


Au  (na)=2|ia   I   £   f~G   (na-pa ,  qa  ,0)  -  G   (na-pa,  qa,  a)~j  (54) 

p=0  q=-°° 

changing  variables  by  setting  u  =  n-p  changes  Equation  (54)  into 

2  /n     °°\  +"  r  -i 

Au2(na)  =  2pa   J  E  +   E  j  E   ^G22  (ua  ,qa,0)  -  G22(ua,qa,  a)J 

\u=l   u=l/  q=-=° 


+    4^2   Z  QG22(0,qa,0)  -  G22(0,qa,a)] 
q=l 

+  cTnT'2^2  G22(°'°'a)'  (55) 

where  C  (n)  is  the  force  constant  defined  in  Equation  (7)  for  a  point 
force  acting  on  a  point  at  a  distance  na  from  the  origin.   After  computa- 
tion of  the  sums  as  presented  in  Equation  (26) ,  Au  (na)  becomes 


n      2 
Au„(na)  = 


2V      C„(n)    2rr 


2 

2.866  -  2.368V    _    .  u   +  1 

-  L    m 7i — 


1  -  v 

u=n+l 


(56) 


For  each  value  of  n,  a  direct  comparison  can  be  made  between 
the  values  of  Au  (na)  from  Equations  (49)  to  (53)  on  one  hand,  and 
Equation  (56)  on  the  other  hand.   Since  both  ought  to  be  identical, 
a  value  of  C„(n)  will  be  obtained  for  each  value  of  n.   Table  2  lists 
the  different  values  of  Au  (n)  and  C  (n)  for  n  =  0  to  5  and  for  n  =  °°. 
The  values  of  C  (n)  are  very  nearly  constant,  and  it  seems  reasonable 
to  expect  a  value  of  C  (0)  very  close  to  C  (°°)  .   This  extrapolation 
permits  us  to  evaluate  Au  (0) . 


TABLE  2.   Variation  of  Force  Constant  C  with  the 
Atomic  Positions  in  a  Screw  Dislocation 
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Aur 


-(n) 


2^ia 
C„ 


Au, 


'(X1)      C2(n)/\iA 


2885 
3988 
4343 
4511 
4607 
4670 

4957 


.9033 

.3964 

.9387 

.3965 

.9554 

.3966 

.9651 

.3966 

.9713 

.3966 

.0000 

.3966 

Replacing  C  (n)  by  C  (°°)  in  Equation  (56)  for  any  values  of  n 


and  noticing  that  for  large  n 


.   .        2iia      a  2.866  -  2.368V 
AU2(na)  =  3  =  cT(nT  +  2^  T^ 


(57) 


a  very  simple  expression  for  the  relative  displacement  across  the 
slip  plane  for  n  ^  0  is  found: 


2 

a    /        ^  a    T-  n      U   +  1 

AU2(na)  =  3  ~  2^  E  ^  2~ 

n+1     u 


(58) 


In  particular  for  n  =  0,  Au  (na)  =  .794a.   This  value  fits  very  well  on 
the  interpolated  curve  for  Au  (x  ) ,  as  shown  in  Figure  8.   For  the 
atomic  points  corresponding  to  x  =  -na ,  with  n  >0  ,  a  direct  transfor- 
mation of  Equation  (48)  leads  to 


AU2(-na)  =  2^  Z  U  2 


(59) 


31 


It  is  striking  to   notice  that  the  symmetry  of  the  screw  dislocation 
displacement  field  is  preserved  at  the  atomic  points.   The  dislocation 
line,  in  the  continuum  sense,  lies  exactly  at  x  =  -a/2.   This  result 
could  have  been  guessed  earlier  by  simply  Considering  an  oriented  path 
around  each  loop.   By  adding  the  loops  together,  the  only  remaining  part 
of  the  path  would  be  a  straight  line  at  x  = -a/2.   This  method  can  be 
generalized  for  determining  dislocation  lines  in  more  complicated  cases. 
This  symmetry  does  not  appear  in  Equations  (49)  to  (53)  because  of  the 
divergence  of  these  expressions  at  x  =  0. 

A  mapping  of  the  atomic  displacements  in  atomic  planes  imme- 
diately above  and  below  the  slip  plane  is  shown  in  Figure  9. 

The  width  of  the  screw  dislocation  is  defined  to  be  the  region 
in  which  the  relative  displacement  is  comprised  between  a/4  and  3a/4. 
Since  the  region  where  the  relative  displacement  is  equal  to  3a/4 
cannot  be  known  exactly,  the  value  2w  =  .73a  for  the  width  can  be 
obtained  by  rough  measurements  on  Figures  7  and  8. 

The  relative  displacement  described  above  can  be  compared  to 
the  results  obtained  by  Peierls  [20,21].   Fitting  the  expression  for 
relative  displacement  obtained  from  the  Peierls   model  to  our  results 
leads  to  the  expression 

-1  xi  +  d 

Au  (x  )  =  |  +  -  tan   — ,  (60) 

2   1    2   TT         w    ' 

where  w  is  the  half-width  and  d  a  translation  parameter  fixing  the 
center  of  symmetry  of  Au  .   For  large  x  ,  this  expression  can  be 
expanded  in  consecutive  powers  of  1/x  , 
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a  w    a   wd  ,  _  .. . 

AU2(X1}  =  a  -^3T+  n+T  -  •••  •  (61) 

1       Xl 

The  same  series  expansion  holds  for  Equations  (49)  to  (53)  and  leads  to 

2 

AW  =  a"2^x-+  4i^+  •••  •  (62) 

1       Xl 

By  comparison  of  Equations  (61)  and  (62) ,  the  half-width  and  the 

translation  parameter  are 

w  =  |  ;    d  =  |  .  (63) 

These  values  correspond  to  the  Peierls'  model  when  the  origin  is  taken 
at  x  =  -a/2.   Although  the  width  of  the  dislocation  in  our  model  is 
slightly  smaller  than  in  Peierls   model,  the  two  expressions  for  the 
relative  displacement  across  the  slip  plane  are  exactly  the  same  for 
large  values  of  x  ,  as  shown  in  Figure  8. 

Finally,  following  Eshelby's  suggestion  [1  ],  a  distribution 
function  for  infinitesimal  dislocations  in  the  glide  plane  can  be  defined. 
Instead  of  being  the  result  of  a  singularity  concentrated  on  the  x   axis, 
the  straight  screw  dislocation  is  considered  to  be  composed  of  a  contin- 
uous distribution  of  infinitesimal  dislocations.   This  distribution 
function  is,  in  fact,  the  component  a       of  the  dislocation  density  tensor 
as  defined  by  Kroner  [2],   It  is  to  be  found  equal  to 

d(Au  ) 
*22(X1>  =  -dx7~  •  (64) 

Differentiating  Equations  (49)  to  (53)  with  respect  to  x  leads  to 
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2  2 

*22   =    2^ 2— 


(65) 


1    ,      n  1 


22v    V    '         2TT         ,22' 
x  (x  +a  ) 


(66) 


and 


2    i      ^  a 


»no(xJ 


2  2 

1     Xl    -   a 

2*      2      2   2 
L_xi    (xi+a   )   _ 


4  2    2         4 

x     -  6a   x    + a 


22      1'         120tt 


jl_ 

2        2.4  4 

L.         ^l^     }  Xl-» 


(67) 


(68) 


It  can  be   easily  verified  that 


+  00 

I     ^2<V    dXl  = 


(69) 


A  plot  of  the  successive  approximations  is  made  in  Figure  10, 
emphasizing  the  values  of  a      (x   )  at  the  atomic  positions.   Symmetry 
of  the  distribution  function  at  these  points  with  respect  to  x  =  -a/2 
is  evident. 

In  summary,  we  shall  emphasize  that  this  model  is  in  perfect 
agreement  with  the  previous  techniques  employed  for  obtaining  the  dis- 
placement field  of  a  straight  screw  dislocation.   An  improvement  has 
been  made  in  the  present  case,  obtaining  a  simple  analytic  expression 
for  the  atomic  displacements  around  the  defect  without  any  exceptions. 
Such  atomic  displacements  can  be  obtained  by  lattice  dynamics  computa- 
tion, but  these  have  the  disadvantages  of  being  difficult  to  use  and 
being  an  entire  numerical  method  [13]. 
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Self-Energy  of  the  Screw  Dislocation 

As  it  has  been  explained  in  the  previous  chapter,  such  an  array 
of  point  forces  can  simulate  a  system  composed  of  two  infinitely  long 
parallel  screw  dislocations  of  opposite  signs  if  the  length  La  is  taken 
much  larger  than  the  separation  distance,  or  width  of  the  array,  Ra. 
In  this  case,  the  energy  of  the  system  per  unit  length  of  screw  disloca- 
tion becomes  W/La  when  L  is  large. 

2 


_W_  _  R+l    2 
La    4tt  ^ 


4tT^+A 

2 


\&  R 

2tt    r'  ' 


(70) 


where  the  constants  A  and  r'  are  defined  by 


A  = 


2.8545  -  2.3676V 
1  -  V 


(71) 


and 


&n—  =   1.179 

ro 


(72) 


Taking  the  special  case  of  V  =  1/3  leads  to 


_W  _  R+l   2 

La    4tt  ^ 


4tt  ££  +  3.0979 
C2 


lia 


&r. 


R 


3076 


(73) 


As  already  mentioned,  this  energy  is  higher  than  the  energy  of 
the  system  composed  of  two  parallel  screw  dislocations  of  opposite  sign, 
because  of  the  nonrealistic  strain  energy  stored  in  the  region  between 
the  planes  of  forces.   The  region  of  the  continuum  where  the  strain  is 
larger  than  1/2  is  shown  in  Figure  11  and  has  the  following  boundaries: 


-  <  x  <  - 

2  ~  3   2 


(74) 


•.35a  <  x  <  Ra  +  .35a  . 
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The  points  x  =  -.35a  and  x  =  Ra  +  .35a  are  those  where  Au  (x  )  =  a/2. 
Outside  these  limits,  Au  (x  )  is  smaller  than  a/2.   The  correction 
energy  which  has  to  be  subtracted  from  W/La  is  composed  of  the  strain 
energy  per  unit  length  of  screw  dislocation  due  to  the  displacement 
field  u(x  ,x  ),  minus  the  strain  energy  per  unit  length  of  screw  dis- 
location due  to  the  displacement  field  measured  relative  to  the  final 
atomic  positions  inside  the  slab,  v(x  ,x  )  described  as  follows: 

0 


V1(X1*X3) 


/v2(xi,x3)  =  u2(x1(x3)  -  x3 


(75) 


v3(xlFx3)  =  0 


The  stress  fields  associated  with  these  two  displacement  fields 


^2 

ai2 =  v.  ^ 


du 
2 

723  "    *  o^ 


(76) 


and 


^2  ^2 

Ti2  -  *  sq =  »-sq 


SV2  /3U2         \ 


(77) 


The  correction  energy  is  the  following  integral  computed  over  the  volume 
mentioned  above. 


E  =|  J"   (a..  eU.  -r..  eM  dV 
c   2  J   \  ij   ij     ij   ij/ 


(78) 
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Using  Green's  theorem,  this  volume  integral  can  be  transformed  into 


a  surface  integral , 
1 


i  (o. .  u.  -  T  .  v  )  n   dS , 


(79) 


where  u  is  the  normal  to  the  surface  S.   Using  Equations  (75)  and  (76) 
leads  to 


2  J 

S 


'  2 


ov. 


l2  ^q  ~   v2  o^j  !1 


n.  +  u 


2  oxr 


°X3/   3 


dS    (80) 


u-  r 


du 


sL3dxl  x 


3  ox     3     2 1   3 


(81) 


where  the  surface,  S,  is  composed  of  the  areas 

S  :   x  =  ±  -,    -35a  <  x  <  Ra  +  .35a; 


S2:   xx  =  -.35a, 


S  x„  <  -: 

3    2' 


S  :   x  =  Ra  +  . 35a, 


^  xn  <  - 
3    2 


Taking  the  symmetry  with  respect  to  x  =  Ra  '2  into  account,  and  the  fact 

that  the  integrand  is  not  dependent  on  x  ,  lead  us  to  the  final  formal 

expression  for  E  /La: 
c 

E         a/2     du9 

u  =  -  »I  „  x3^r  ("-35a>   V  dX3 

-  a/2       1 


Ra/2 
-  .35a 


a    2      a     a      ,    a. 

2  ^  (X1'2}  "  2  +  U2(Xr2} 


dx. 


(82) 


The  analytical  expressions  for  u  ,  du  /dx  ,  du  /dx  are  obtained 


from  Equation  (31)  and  its  derivatives: 
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VW 


R 

Z  2n 


<  ~,2  /    ^2 

(x  -pa)   +  (x3+2> 


p=0    (xi~pa)   +  (x3  --) 


(83) 


h   3x~  (X1'X3)  =  ^  Zn  X3 
1  p=0 


and 


x  -  pa 


^    ,    *n2    ,  x2   /     ax2 

(Xl-pa)   +  (x3+-)     (Xl-pa)  +(x3-^)J 


9U2  a   R 

ox„    13    2tt 


p=0 


a 
X3  "2 


t  x2    /    ax2    ,      n2    ,     ax  2 

(x^pa)   +  (x3+-)     (xrpa)   +  (x3--) 


(84) 


(85) 


The  mathematical  problem  of  integrating  a  function  over  a  region  where 
the  integrand  has  singular  points  for  certain  values  of  x   is  removed 
by  the  same  argument  as  the  one  used  in  the  previous  section. 
After  integration,  Equation  (82)  becomes 


pa 


2   R 


La 


p=0 


(p  +  .  35)  2m   (P+'35)  t   1  +  2  tan"1  (p+.  35) 


(P+.35) 


-  V* 


(•-I)- 


(86) 


The  summation  over   p  is  computed  in  the  same  way  as  before  using  Equa- 
tion (26) ,  which  gives  finally  for  the  correction  energy 


E 
c 

La 


[13. 


R 


\» 


rr    r 


R 


(87) 


and 


0n  rJJ  =  .3129. 


(88) 


The  total  energy  of  the  system  composed  of  two  antiparallel 
screw  dislocations  is 
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_T  _   _W_ c 

La  _  La        La 


(89) 


T  _   R+l        2 

La  "     4tt    ^ 


4tt  ^+  A-2tt 
-  2 


2       R-rl  2 

(ja     „_         0        [ia 


2tt 


&■ 


<'o>2      2 


(90) 


From  the  usual    continuum  theory   of   dislocations,    this    energy  has   to  be 
equal    to   the   sum  of   the   self -energies   of   both    screw   dislocations,    minus 
their   interaction  energy 


La 


ua 
4rr 


\» 


2m  -  . 


2rr        R 


(91) 


The  requirement  that  our  expression  (90)  has  to  be  identical  in  form 
to  Equation  (91) ,  will  force  us  to  choose  the  still  unknown  parameter 
C_  such  that  the  term  divergent  with  R  vanishes.   Thus, 


with 


if 


E      2 

La    2tt    r 


2m  —  =  1.3367, 

ro 


|ia  _  1 


2.8545  -  2.3676V 
4tt(1-V) 


For  V  =  1/3 


T  _  _ya_ 
La  ~  2tt 


2r. 


Ra 


2631a 


(92) 


(93) 


(94) 


(95) 


if 


C  =  3.945(ia 


(96) 
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From  Equation  (91) ,  the  self-energy  of  a  pure  single  screw  dislocation 
can  be  written 

E„     2 

r~  =   *T=-fa  .   .  (97) 

La    4tt    .  263a 

Two  important  remarks  can  be  made  here.   First,  the  value  found 
for  C   is  very  close  to  those  found  from  the  displacement  field  computa- 
tions.  This  proves  the  consistency  of  the  correction  energy  with  the 
displacement  field  included  from  the  array  of  forces.   Secondly,  it  is 
found  that  the  core  parameter  in  Equation  (97)  is  independent  of  V, 
Poisson' s  ratio.   This  is  in  complete  agreement  with  Peierls'  result, 
which  gives  a  value  of  r   equal  to  a/e,  that  is,  equal  to  .37a,  where 
e  is  the  naperian  base  of  logarithms. 

So  both  models  give  nearly  identical  results,  but  with  a  slightly 
different  r  . 

For  a  direct  comparison  with  Volterra' s  dislocation  model,  the 
core  radius  r   in  Volterra' s  model  has  to  be  reinterpreted  and  cannot  be 
considered  anymore  as  a  cut-off  core  radius  where  Hooke' s  law  does  not 
apply.   It  is  rather  a  constant  containing  all  the  constant  terms  aris- 
ing in  the  computation  of  the  core  energy.   The  Volterra  cut-off  radius 
could  be  evaluated  in  a  better  way  by  defining  the  region  where  Hooke' s 
law  does  not  apply,  e.g.,  where  the  strain  is  larger  than  .10.   From 
Figure  8,  such  a  cut-off  radius  can  be  approximated  as  being  about 
rQ  =  1.5a. 

So,  the  technique  of  simulating  dislocations  by  a  point  force 
array  seems  to  be  very  successful  in  describing  the  principal  features 
of  the  defect,  even  though  a  complete  accuracy  in  the  computation  of  the 
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atomic  displacements  cannot  be  reached  because  of  the  elastic  and 
isotropic  approximation. 

Single  and  Double  Kinks  in 
a  Screw  Dislocation 

A  single  or  a  double  kink  can  be  simulated  by  simply  adding  to 

the  array  of  forces  an  extra  row  of  shear  loops,  parallel  to  the  screw 

dislocation  line,  and  a  semi-infinite  or  finite  extent,  respectively. 

A  representation  of  the  modified  array  is  sketched  in  Figures  12  and  13. 

The  displacement  field  and  the  energy  of  these  defects  are  handled  in 

the  same  way  as  for  the  straight  screw  dislocation. 

A.   Displacement  field 

The  displacement  field  of  these  defects  is  obtained  by  adding 
to  the  displacement  field  of  the  pure  screw  dislocation,  the  displace- 
ment field  resulting  from  the  extra  row  of  forces.   The  latter,  u' ,  has 
the  following  expressions  for,  respectively,  a  single  kink  and  a  double 
kink  of  length  2Na. 

CO 

q=0 

2 
"   G2i(Va'X2-qa'X3+lQ    "   ^   CG2i(Xl+a'VX3-f} 

2 

-  G2i(va'Vx3+ti]  -  ^Qwv^v^-f* 

+    G3i(Va'X2'X3  +  l}3    •  (98) 
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Figure  12.   Array  of  Forces  for  a  Single  Kink 
in  a  Screw  Dislocation 
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Figure  13.   Array  of  Forces  for  a  Double  Kink 
,  in  a  Screw  Dislocation 
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2      N 


U1(DK)(X1'X2'X3)    =    **         E        C^i^l^'V*18'^-^ 


q=-N 

2 

-  G2i(xi+a-vqa'x3+l9  ■JrCG2i(v,'yM'x3"3) 

-  G2i(Xl+a'VNa'X3  +  l)    +    Si^l^'V^'^"^ 

2 

-  G21(Va'VHa'X3  +  l)J  +  Jf-CG31(Va'VHa'JI3-|) 

+  G3i(xi+a,x2-Na,x3+J)  -  Gg.  (Ya'VNa'X3  "  |} 

-  Si'V^V^s  +  M  •  (99) 

Computations  of  u'  and  u'  for  both  kinds  show  that  these  displacements 
are  very  small,  even  of  the  region  of  high  distortion  in  the  x  direc- 
tion.  They  reach  a  magnitude  of  a  few  thousandths  of  an  atomic  distance. 
Therefore,  we  shall  concentrate  our  attention  on  the  u'  component  of  u' , 
and,  more  specifically,  on  the  atomic  displacements  in  the  planes  just 
above  and  below  the  slip  plane  (x  =  ±  a/2) ,  since  this  is  the  region  of 
highest  distortion. 

Because  of  the  existence  of  singular  points  at  the  points  of 
application  of  the  point  forces,  several  special  cases  will  be  consid- 
ered. As  a  first  step,  we  shall  restrict  our  range  of  interest  by  notic- 
ing the  various  symmetries  in  the  expression  of  u' .   It  shows  an  odd 
symmetry  with  respect  to  x  =  0 ,  and  an  even  symmetry  with  respect  to 
x  =  -a  for  both  cases,  single  and  double  kinks.   It  shows  an  even 
symmetry  with  respect  to  x  =  0  in  the  special  case  of  the  double  kink. 

(1)  Displacement  field  of  the  single  kink 

The  single  kink  will  be  the  first  case  considered.   For  most 

of  the  values  taken  by  x_  and  x„ ,  u '  . _  .  has  the  following  expression: 

1      Z        <s(SK) 
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'  (  ai    _      a 

U2(SK)U1'X2'2;    -    8^ 


Z      f(x2-qa)    -   -  f(x2) 
Lq=0 


2 
a   x„ 
_a 1_       2 

32tt  1-v      r  o        .         r,13/2 


T  .2        2        2T 

UCx  +a)     +  x    +  a  _J 


(100) 


with 


f(X2-qa)    =   2(1^0 


/  2  2  /  2  2  2 

V(x  +a)      +    (x  -qa)  V(x  +a)      +    (x  -qa)      +   a  _ 


2(l-v) 


(x2-qa)' 


(x2-qa) 


f  2  2~|3/2       r  2  2        2i3/2 

_L(x  +a)      +    (x   -qa)  J  (Jx  +a)      +    (x  -qa)      +    a  J 


(101) 

On  atomic  positions,    different    expressions    apply  because  of   the   singular 
points   situated  on  x     =   -a: 

(a)      x     =   na      and     x     /  -a 


U2(SK)(Xl'na'l)    =^ 


»  n 


E     +      E        f  (qa)    -   »   f  (na) 
q=l        q=0v' 


32rr(l-v) 


C2  2    2  2  I 

x„  +  a)      +   n  a     +    a  J 


3/2 


(102) 


(b)      x     =    -na      for   any   x 
«  1 


^(SK)^!'""3'^    =   ^ 


E     -      E        f(qa)    +   £  f(na) 
L\q=l        q=-  ' 


32tt(1-v) 


f,  ,2  2    2  2~] 

MX  +a)      +    n  a     +   a  -J 


3/2 


(103) 
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(c)      x     =   na     with      n  >  0,    and   x     =    -a 


,2<siOC",,l'na'3) 


Hi-iX-&-7=)^ 


vq  +j 


■V)  3/2 

(q   +D  J 


2i  ! 


n+1 


5~r  2(l-v)  3/2 

(n  +1)        ~ 


Jxf+1 


|ja    _   a      3-4V  . 

+  ~cj   16rf    1-V 


104) 


(d)      x     =   0      and     x     =    -a 


q=l 


2 (1-V)  3/2 

(q  +D        -' 


pa  a       /3-4v\ 

2C~       32rf    \  1-V  J    ' 


(105) 


In  all  these  cases,  f  has  the  same  form  as  in  Equation  (101). 
Q.     has  the  same  meaning  as  before.   Since  it  has  been  noticed  that  the 
force  constant  is  very  nearly  constant  for  all  rows  of  point  forces,  we 
shall  give  it  the  value  found  from  Equation  (57).   The  atomic  displace- 
ments are  tabulated  in  Table  3  for  the  following  values  of  x   and  x  : 
x  =  -2a,  -a,  0,  a 

x  =  -5a,  -4a,  ...,  0,  ...,  4a,  5a  . 

To  these  values  have  to  be  added  the  corresponding  atomic  dis- 
placements due  to  the  straight  screw  dislocations.   A  mapping  of  the 
atomic  configuration  above  and  below  the  slip  plane  is  shown  in 
Figure  14.   This  sketch  shows  a  very  high  distortion  around  the  point 


(x  =  -a,  x  =  0) ,  but  after  a  few  atomic  distances  from  this  point, 
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the  atomic  configuration  shows  very  little  difference  from  the  configu- 
ration of  the  pure  screw  dislocation.   Even  if  the  actual  width  of  the 
kink  cannot  be  expressed  analytically,  one  can  conclude  that  the 
defect  is  very  localized. 


TABLE  3.   Atomic  Displacements  for  a  Single  Kink  in  a 
Screw  Dislocation 


x2/a 

u' 
2(SK) 

/a 

x  =  -2a 

xi=  -a 

x  =  0 
1 

xi=  a 

-5 

.0019 

.0020 

.0019 

.0016 

-4 

.0028 

.0030 

.0028 

.0023 

-3 

.0046 

.0052 

.0046 

.0034 

-2 

.0086 

.0111 

.0086 

.0052 

-1 

.0173 

.0341 

.0173 

.0073 

0 

.0347 

.1473 

.0347 

.0102 

1 

.0374 

.2605 

.0374 

.0105 

2 

.0461 

.2834 

.0461 

.0126 

3 

.0501 

.2893 

.0501 

.0163 

4 

.0519 

.2915 

.0519 

.0154 

5 

.0529 

.2926 

.0529 

.0161 

(2)   Displacement  field  of  the  double  kink 


We  shall  consider  here  a  double  kink  of  length  2Na.  The  compu- 
tations have  the  same  form  for  a  kink  of  length  (2N+l)a.  The  displace- 
ment field  at  a  general  point  in  the  plane  x  =  a/2  has  the  following 


expression: 


,2(DK)(X1,X2,2) 


2 

a 

8Tr 


N 
E 

_q=-N 


32rr(l-v) 


f (x2-qa)    -  -  f(x2-Na)    -  i  f(x2+Na) 


x  +Na 


_Gv 


,3/2 


2  2      2~Y  f 

a)    +(x  -Na)    +a  J  [_(x  +a)    +(xn+Na)    +a 


2  .2      2"| 


3/2 


(106) 
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with 


f(X2"qa)  =  2(1^) 


I  o         2      '      2 

j/(x  +a)  +  (x  -qa)     */(x  +a)  +  (x2~qa) 


2   2 
+  a 


2(1-V) 


(x2-qa) 


(x2-qa) 


r  2  2i3/2      T    2      2  2~f/2 

Jjx  +a)  +(x2-qa)  J      [jx^  +  a^+U^qa)  +a  J 


(107) 


The  points  where  this  expression  do  not  apply  are  the  points  of 

application  of  the  point  forces  on  the  row  x  =  -a.   Two  separate 

cases  are  considered. 

(a)      x     =   -a      and      x     =   ±  na      for      n  <  N 
1  & 


^(DK)*-*!*18'^    =    8t 


E     +      E     ) f(q)    -  ~  f(N+n) 


32tt(1-V) 


\'+n        N-n 
E     + 
q=l       q 

N-n  N+n 


JjN-n)2+lJ  QN+n)2+lJ 


2 

pa  a 

C2    -16 


|  f  (N-n)J 


(108) 


(b)      x     =    -a      and     x     =   ±  Na 


ui(K)<-a'iNa'2)    =    8^ 


2N 
E 
Lq=l 


E      f (q)    -   2    f(2N) 


16tt(1-V)  2        3/2 

(4N  +1) 


pa  a         /3-4\) 


2C  32rr      I  1-v  /  ' 


(109) 


with,    in  both   cases, 


f(q)    -    2[    i L 


-/q2+: 


2  (1-V)  3/2 

(q   +D 


(110) 
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The  sizes  of  double  kinks  have  been  considered,  corresponding 

to  the  values  of  N  equal  to  1,  2,  and  3.   Due  to  the  symmetry  already 

mentioned,  only  the  values  of  u'      for  the  following  values  of  x  and 

2 ( DK )  1 

x  have  been  computed: 

x  -  -a,  0,  a,  2a 

x9  =  0,  a,  2a,  3a,  4a,  5a  . 

These  values  of  u'      are  tabulated  in  Tables  4,  5,  and  6, 
2d  (  DK  ) 

corresponding,  respectively,  to  N  =  1,  2,  and  3. 


TABLE  4.   Atomic  Displacements  for  a  Double  Kink  of  Length  2a 
in  a  Screw  Dislocation 


x2/a 

U2(DK)/a 

x  =  -2a 

1 

x  =  -a 
1 

xi  =  0 

*1  =  a 

x  =  2a 

0 

.0201 

.2374 

.0201 

.0031 

.0009 

1 

.0188 

.1363 

.  188 

.0037 

.0011 

2 

.0127 

.0288 

.0127 

.0038 

.0014 

3 

.0058 

.0081 

.0058 

.0028 

.0013 

4 

.0028 

.0033 

.0028 

.0018 

.0011 

5 

.0015 

.0012 

.0015 

.0011 

.0008 

TABLE  5.   Atomic  Displacements  for  a  Double  Kink  of  Length  4a 
in  a  Screw  Dislocation 


U2(DK)/a 

X2/a 

x  =  -2a 

x  =  -a 

1 

xi=0 

Xl  =  3 

x  =  2a 

0 

.0375 

.2781 

.0375 

.0074 

.0022 

1 

.0328 

.2661 

.0328 

.0080 

.0023 

2 

.0246 

.  1444 

.0246 

.0066 

.0024 

3 

.0155 

.0321 

.0155 

.0057 

.0024 

4 

.0073 

.0097 

.0073 

.0040 

.0020 

5 

.0037 

.0042 

.0037 

.0005 

.0006 
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TABLE  6.   Atomic  Displacements  for  a  Double  Kink  of  Length  6a 
in  a  Screw  Dislocation 


x2/a 

u'     /a 
2(DK) 

x  =  -2a 

x1=-a 

x  =  0 

1 

xi=a 

x  =  2a 

0 

.0455 

.2872 

.0455 

.0109 

.0036 

1 

.0429 

.2861 

.0433 

.0103 

.0035 

2 

.0356 

.2694 

.0356 

.0088 

.0033 

3 

.0261 

.1461 

.0261 

.0077 

.0031 

4 

.0164 

.0330 

.0164 

.0064 

.0029 

5 

.0089 

.0103 

.0089 

.0044 

.0024 

A  mapping  of  the  atomic  arrangement  is  attempted  in  Figures 
15,  16,  and  17,  corresponding  to  the  cases  N  =  1,  2,  and  3,  respectively. 
The  same  remarks  can  be  made  about  the  double  kink  configuration  as  has 
been  made  for  the  single  kink.   Once  again,  the  kink  seems  to  be 
a  very  localized  defect  for  a  structure  like  a  simple  cubic  crystal. 
In  the  case  of  N  =  3 ,  that  is,  where  the  length  of  the  double  kink  is 
6a,  each  single  kink  part  of  the  whole  kink  seems  to  behave  like  a  pure 
single  kink,  which  is  to  be  expected  when  the  double  kink  grows  in  size. 
This  argument  will  be  used  to  extrapolate  the  energy  of  a  single  kink 
from  the  energy  of  a  very  long  double  kink. 


B.   Energy  of  a  double  and  a  single  kink 

The  energy  of  a  double  kink  is  defined  as  being  the  difference 
between  the  energy  of  the  modified  array  of  forces  and  the  energy  of  the 
rectangular  dislocation  loop.   As  has  been  developed  for  the  screw 
dislocation,  a  correction  energy  term  will  be  introduced  to  take  into 
account  the  actual  strain  field  across  the  slip  plane. 
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In  a  first  step,  the  difference  between  the  energies  of  the 
systems  of  forces  has  to  be  computed.   It  is  simply  the  self -energy  of 
the  extra  double  row  of  forces  of  length  2Na,  and  the  interaction  energy 
between  this  double  row  of  forces  and  the  dislocation  loop  array. 


W.JLL.  (4N-1>|  A-G2a(0>0fa) 


2  4 

\x   a 


-+G33(0,0^ 


2  4 


+  —~-   (jG22(0,2Na,0)  -G22(0,2Na,a)  -  G33<0,2Na,0)  -  G33(0 ,  2Na,  a)) 


2  4  2  4  2N~1        r  1 

-  |*  a    G93(0,2Na,a)  +  2\i   a    2    (2N-q)  |^G  2(0,qa,0)  -  G  2(0fqa,a)J 

q=i 

2  4  2N_1  2  4  R+1  f  La 

■2|iB   2    G23(0,qa,a)  +  (i  a   2  <  G22(pa,— -  Na,0) 

q=l     "  p=l  L 


-  G22(pa,-^-Na,a)  +  G22  (pa  ,-— +  Na,0)  -G^Cpa,-^  Na,a) 


+  G   (pa.-^-Na.O)  +  G33(pa,^  -  Na,  a)  ~G33(Pa'-^+  Na,0) 


La  La  L/2+N-1 

-  G33(pa,— +  Na,a)  -  2G23(pa,-2^-+ Na,a)  -22      G23(pa,qa,a) 


+  4N  j^G22(pa,0,0)  -  G22(pa,0,a)T) 


L/2-N-1     L/2+N-1 
8N  2     +22     (L+2N-l-2q) 

q=l 

"  "2 


q~-N 


Q}22(pa,qa,0)  -  G22(pa,qa,  a)J 
(111) 


We  shall  consider  cases  for  which  the  length  of  the  double  kink 
is  much  smaller  than  the  length  of  the  dislocation  line.   Under  this 
condition,  t',       expression  of  W  is  greatly  simplified. 
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(2N-1 
V*       1    3-4v     ,    1       v 
el     8tt  2(1-v+4tt     L 
2 

q=l 


,  ,l._i 


«/q  +1 


2(1"V)     (q2+l)3/2J 


_1_ 

8" 


q=l 


3-4V  .    q  +1  IIS 

T^7^~2-+  1^1  "2 


q  +i 


l 

—  + 
q 


Vq2+1 


,    2         3/2 

(q  +D       j 


pa    i     pa         1        1  1 

2C    +2C^+32rr  1-v     32rr 


7-8v 


4N 


1-v 


-s/4N  +1 


1-v  2        3/2 

(4N   +1) 


2N-1 

2^     S 
q=l 


1    - 


/  2 
Vq   +1 


2(1-V)     ,    2        3/2 
(q  +1) 


J 


(112) 


The  force   constant   C     has   been  determined   from  previous    computations 
for   the   screw  dislocation.      Assuming   that   the   slight  modification  of 
the  array  does  not  have  any   influence  on  the  value  of   C2,    it   is   deter- 
mined  such   that 


|ia         A    _    1 
S   +   4n  =   2    ' 


(113) 


with 


1      3-4V     _1_  3-4V  /l 

16rr    1-V  +  8tt    1-v         ,      q 

q=i  l 


1  N 

\          1     3-4V      "     „     q2+l 
+    8TT     1-V       =    *»  -T- 

'                    q=i        q 

</q2+i  J 

.      1       1        y 

+  4rr  1-v     '- 


q=l   q   +1 


(114) 


So,    W/pa     becomes,    after   transformation, 


p,a  I  2N 


/  2 

vq  +1 


2(1-V)         2         3/2 
(q  +1) 
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\ia         |ia  1        1  1 

2C^  +    2C^  +    327f  1-V        327f 


7-8v    1 
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Jsi 
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4N 


3/2 


4N  +1 


(4N  +1) 


2N-1 


2rr 


q=l 


^ .    2(1-V)    -2    3/2 

</q  +1  /  (q  +1) 


(115) 


If  Na  is  sufficiently  large,  but  still  smaller  than  the  length 
of  the  dislocation  line,  the  corresponding  double  kink  will  behave  like 
two  separate  single  kinks.   Each  will  have  an  energy  equal  to  half  the 
energy  of  the  whole  kink,  that  is, 


W 


\ia 


N      -  ^  P-3  .  1  1  1  y 

2C     +  2C0   +   32tt  l-v"2rr        , 
2  3  q=l 


1    - 


/2~ 


2  (1-V)         „        3/2 
(q  +1) 


(116) 


The  difference  between  twice  the  energy  of  a  single  kink  and  the  energy 
of  a  double  kink  represents  the  interaction  energy  between  the  kinks. 
Two  remarks  are  necessary  here.   First,  a  term  proportional  to  the 
length  of  the  kink  is  included  in  the  final  expression,  and  is  expected 
to  cancel  out  with  a  similar  term  in  the  correction  energy.   Secondly, 
the  force  constant  C„  appears.   This  has  to  be  considered  an  unknown 
parameter,  since  there  is  no  physical  condition  which  can  be  applied  to 
evaluate  it.   An  approximation  could  be  made  by  setting  C   equal  to  C  , 
but  there  seems  to  be  no  particular  justification  for  such  an  assumption. 

The  second  step  needed  to  obtain  the  final  expression  for  the 
energy  of  a  double  kink  is  the  computation  of  the  difference  between 
the  correction  energy  for  the  pure  screw  dislocation  and  the  correction 
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energy  for  the  kinked  screw  dislocation.   Since  the  former  has  already- 
been  computed,  we  shall  focus  our  attention  mainly  on  the  latter. 

The  region  that  suffers  a  strain  larger  than  one-half  is  bounded 
by  the  planes  x  =  ±  a/2,  the  surface  parallel  to  x  where  the  relative 
displacement  is  equal  to  a/2  and  the  plane  x  =  Ra/2  (see  Figure  18). 
The  region  situated  between  x  =  Ra/2  and  x  =  Ra  will  not  be  con- 
sidered, since  its  deformation  is  the  same  as  for  the  straight  screw 
dislocation. 

We  have  already  seen  that  the  displacement  field  caused  by  the 
array  of  forces  is 


'"l^i-W 


swv +  ^wv 


(117) 


V.U3(VW  • 


By  taking  into  account  the  fact  that  after  relaxation,  the  atomic 
positions  need  to  be  referred  to  their  closest  neighbors  (see  Figure  19) , 
we  are  led  to  choose  as  actual  displacement  field  across  the  slip  plane 
the  following  expressions: 


^v    (x  ,y  ,x„)  =  u' (x  ,x„,x„) 
lv  l'J2'  3      1   1'  2'  3 


\   V2(Xl'y2'X3)  =  U2(X1'X3)  +  ^l'W  ~  X3 


(118) 


v  v„(x  ,y„,x„)  =  u' (x  ,x„,x„) 
V.  3   l'J2'  3     3   1'  2'  3 


where  v      is   related   to  x      such   that 
J2  2 


x3>0 


x3<0 


2    ' 
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_Rq/2. 


t 


x. 


Figure  18.   Region  of  High  Strain  for  a  Double  Kink 
in  a  Screw  Dislocation 


Figure  19.   Atomic  Relaxation  for  a  Double  Kink 
in  a  Screw  Dislocation 
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The   stress   fields   corresponding   to  u   and  v  have   the   following 
components   a.  .    and  T. .,    respectively, 


dug  /^        ou3 

CT22(X1'X2'X3)    =    (X+2^    o^  +    M^  +    o^ 


33 v    1'    2'    3' 


J2         Su2         3U1 

J12(JVW     ^  [o^+  o¥^+   o^ 


ff13(Xl'X2'X3} 


1        ^3 


V 


2         BU2         3U3 

CT23(xi'x2'x3)  =  li\^q  +  ^q  +  ^q 


(119) 


and 


Tll(xi'y2'X3)    =    'WWV 


WW    =    WWV 


t33(x   ,y    ,x  )   =  a33(x1,x2)x3) 


T12(Xl'y2'X3)   =   ^VW 


i*  ~(x    ,y„,x„)   =  a  „(x   ,x„,x„) 
13v    1,J2'    3  13V    1'    2'    3 


T23(Xl'y2'X3)    =    a23(X!'X2'X3)    "    » 


(120) 


Recalling  Equation    (79)  ,    the  general   form  for  the   correction 
energy   is   the   difference  between  the   strain   energies    in  the  region 
of   interest, 
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E  =  ^  f  a.  .  u.  n.  dS  -  -i  f  T.  .  v.  n.  dS'  .  (121) 

c   2  >-s   ij   i   j      2  <^,   ij   i   j 

The  surfaces  S  and  S'  differ  only  by  the  range  of  integration  over  x 
and  y  ,  respectively.   When  x  is  positive,  y   is  defined  between 
-La/2  +  a/2  and  La/2  +  a/2,  and  when  x   is  negative,  y   is  defined 
between  -La/2  -  a/2  and  La/2  -  a/2.   Under  these  conditions,  the  inte- 
grals involved  with  functions  of  y   can  be  considered  as  integrals  of 
functions  of  x  with  different  limits.   The  various  identities  follow: 

La/2+a/2  La/2 

/      f(y  )  dy  =  J  f(x  )  dx    for   x  >0,      (122) 

-La/2+a/2         *        -La/2 

La/-a/2  La/2 

J  f  (y  )  dy  =  J  f (x  )  dx    for   x  <  0  .     (123) 

-La/2-a/2  -La/2 

So  E   can  be  written  as  a  surface  integral  over  a  function,  depending 
on  a  unique  variable  x  . 

E  =  i  f  (CT. .  u.  -  T. .  v.)  n  dS  .  (124) 

c   2  ^   ij   i    ij   i 

Replacing  T. .  and  v.  by  their  expressions  as  functions  of  a        and  u  , 
ij      i  ij      i 

respectively,  leads  to 

E      = 


\  I  (X3   C21   "l   +    (X3CT22    f    ^3)n2   +  OvV  +  X3CT23]  n3)  dS » 


(125) 


i  V2 

Ec  =   2  J     /0   dX3  /CX3CT21  ni  +    (X3C22  +    ^3)n2]dS 
-a/2  C 

a/2      La/2 

"J    dX3  /    X3G21  (^T'X2'X3)  dX2  dX3  + 
-a/2   J  -La/2   J  Z1       "       Z      J    2        3 
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La/2  Ra/2 

+     /    ,,   dX2     J\      .     GVwt*    "    •*  I  +   I  ff23(xl'Vf}J   dXl    • 

-La/  .s     X  v^Q'' 

(126) 
In  these  integrals,  C  is  the  curve  defined  by 


Ws*  +  "J^i'Va?  "  4  •  (127) 


that  is,  the  "line  of  the  kink,"  and  X  (x  )  corresponds  to  a  point 
(x  ,x  )  on  this  curve. 

The  evaluation  of  these  integrals  is  a  complicated  mathematical 
problem  which  has  to  be  solved  numerically.   Unfortunately,  such  numer- 
ical computations  have  not  been  possible  to  achieve  yet,  mainly  because 
of  the  very  complicated  expressions  for  the  displacement  and  strain 
fields.   However,  further  research  on  this  mathematical  problem  can  be 
carried  out  and  will  lead  to  the  correct  answer  for  the  energy  of 
a  double  kink.   The  final  expression  for  the  energy  would  be  obtained  by 
subtracting  the  difference  between  the  correction  energies  given  by 
Equation  (115).   The  linear  term  is  expected  to  cancel,  so  that  the 
energy  is  a  finite  number.   The  limit  of  this  number  when  N  becomes 
large  would  be  twice  the  energy  of  the  single  kink. 


CHAPTER  5 


EDGE  DISLOCATION  IN  SIMPLE  CUBIC  CRYSTAL  CONSTRUCTED 
FROM  AN  ARRAY  OF  SHEAR  LOOPS 


The  same  procedure  as  followed  in  the  case  of  the  screw  disloca- 
tion will  be  used  in  the  case  of  the  edge  dislocation. 

Displacement  Field 

As  seen  previously  in  the  chapter  concerning  the  whole  rectan- 
gular array  of  point  forces,  the  regions  where  the  dislocation  loop  has 
a  pure  edge  character  are  delimited  by 

\    s 

and 

La 
_  —  -  e 
2    2     2 

where  e   and  e   are  small  compared  to  Ra  and  La,  respectively.   Since 
both  regions  simulate  two  identical  parallel  dislocations  of  opposite 
sense,  we  shall  only  consider  the  first  one.   A  translation  of  the  x 
axis  to  the  first  row  of  forces  x  =  -  La/2  will  simplify  the  expres- 
sions.  The  same  symbols  x  and  x  will  be  kept,  having  now  the  meaning 
of  e  and  e   respectively.   This  part  of  the  array  of  forces  is  repre- 
sented in  Figure  20.   The  only  difference  with  Equation  (21)  is  now  the 

L      L-l 
summation  of  q  will  be   Z  and   E   .   So,  u  (r)  becomes: 

q=0     q=l 

2   R/2     j  /l    L-l\ 
V?>  =  "IT    ^  R/9    K  +  \)   CG2m(VPa' V^'X3  -f>  " 

p=  -r/2    n,q=°     q=iy 
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"  'WV^'V^S  +  fD    "  CG3m(xi"pa'3S2'x3-5) 

+    G3m(XrPa'X2'X3  +  IO    +    (?3a(x1"pa'VLa»JI3-|) 
+    G3m(xrpa'VLa'X3  +  H>- 


(128) 


The  summations  on  p  and  q  will  be  accomplished  as  before,  with 
the  help  of  Euler' s  formula  (Equation  (29)).   The  first  summation 
carried  out  will  be  on  p,  since  the  displacement  field  of  the  edge  dis- 
location ought  to  be  independent  of  x  .   Euler' s  formula  becomes  simply, 


R/2  +=> 

2    f(p)  =  /   f (x)  dx  , 

p=  -  R/2        -  co 


(129) 


where   f   is    a   symbol    for   the  whole   expression   to  be   summed.      The  dis- 
placement  field   takes   the  form: 


WW  -  ° 


\q=0        q=lj    1  (x2-qa)      +    (Xg-g) 


1 

1-V 


.  a. 2 

(x3 +  i) 


(X3    "    2} 


(x2-qa)      +    (x3+-)  (x2-qa)       +    (Xg--) 


(130) 


16rr(l-v) 


V^-^  X2(X3+2) 


2  a,2  2  a   2 

X2   +    (x3'2)  X2  +    (X3+2} 


(131) 


Uo(x, ,X       X    ) 


3V    1'    2'"  3'         16tt(1-V) 


L  L-l 

E     +     E 

q=0         q=l 


(x2-qa)(x3--) 

,  %2  t  a^2 

(x2-qa)       +     (x3--) 
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(x9-qa)(x3  +  -) 


a     1-2V  .    R 

S/r. 
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/4+L 


-2  ^2  4"      1-V 
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a    3-4v 


8rr    1-v 
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a      3-4v 
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+  ^<CX2  +  (x3  +  l)2lr   "  155133 


(X3-2> 


a  2 

(X3+2} 


2         *S2     2   ,     <S2 
X2  +  (X3-2}     X2+(X3  +  2> 


(132) 

Unfortunately,  these  expressions  are  much  less  simple  than  those 
found  for  the  screw  dislocation.   But  some  of  the  essential  features  of 
the  displacement  field  of  the  edge  dislocation  can  be  noticed  already, 
that  is,  the  lack  of  displacement  in  the  x  direction  and  the  fact  that 
u   and  u   are  independent  of  x  . 

Euler's  formula  applied  a  second  time  will  give  us  the  final 

expression  for  the  displacement  field.   The  mathematical  difficulty, 

arising  because  of  the  singular  points  in  the  integrand,  is  overcome  in 

the  same  way  as  for  the  screw  dislocation.   The  same  symbols,  u   and 

u  ,  will  be  used  for  the  successive  terms  added  to  approximate  u   and 
3  2 

u  ,  respectively.   For  u  (x  x  x  )  we  find 
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_a 1_ 

16tt  1-v 


V^S"^ 


X2(X3  +  2) 


2  a.2  T      2  a, 2 

L^*^"^  X2+(X3  +  2)    . 


(133) 


1,  ,  a      3-4V 
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(134) 
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(136) 


with 


0  12  3 
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2  2  2  2  2 
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For   u    (x      x    , x   )    we   find 


/2  2 
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Uo(X,  .Xo>Xo)  = 
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(139) 


u„(x    ,x   ,x  )= 
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with 


0  12  3 

U3   =   U3   +   U3   +   U3   +   U3   +    • * • 


(142) 
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The  constant  in  the  expression  of  u  means  only  that  the  point  of 
non-bending  of  the  lattice  planes  is  set  at  x  =  La/2.   However,  the 
relative  positions  of  the  atoms  with  respect  to  each  other  are  inde- 
pendent of  this  constant,  and  the  stress  ana  strain  fields  will  not 
depend  on  the  constant  terms. 

The  component  u  (x  x  x  )  has  four  singular  points  (x  =  0, 
o       \       Z       o  £ 

x  =  a  with  x  =  ±  a/2)  and  u  (x  ,x  ,x  )  has  two  singular  points 
Z,  o  o   X  t*       «j 

(x  =  0  with  x  =  ±  a/2).   This  is  not  surprising,  since  the  whole 

array  of  forces  is  a  superposition  of  two  arrays  of  forces  having 

2 

magnitude  |ja  /2,  one  starting  at  x  =  0  and  one  starting  at  x  =  a. 

An  asymptotic  expression  for  u   and  u   can  be  computed  by 
considering  x   and  x  large  with  respect  to  the  atomic  distance. 
The  following  expressions  are  found: 


u1(x1,x2,x3)    =   0, 
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2      3  2(l-v)      2  2 
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These   equations   can  be   compared   to   those   arising   from  the 
ordinary   continuum  model    for   an  edge   dislocation  of   the   same   sign   [20]  : 
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(148) 


Both  sets  of  equations  are  identical  when  one  is  aware  that  u 
is  determined  only  up  to  a  constant.   Adding  the  expression  a/8rr(l-V) 
to  Equation  (145)  leads  automatically  to  Equation  (148);  this  physically 
means  a  change  in  the  "cut  plane."   One  can  notice  too  that  these 
asymptotic  expressions  for  u  come  only  from  the  first  approximation  u  . 
This  means  that  for  a  point  situated  at  a  large  distance  from  the  dis- 
location line,  the  discrete  array  of  point  forces  appears  to  be  a  con- 
tinuous distribution  of  force  on  the  two  planes  x  =  ±   a/2. 

A  mapping  of  the  atomic  displacements,  except  for  the  singular 
points,  is  shown  in  Figure  21.   Obviously,  this  model  shows  a  strong 
dissymmetry  with  respect  to  the  extra  half  plane  (x  =  .75a,  with 
x  <  0)  in  the  region  of  the  core. 

The  relative  displacement  across  the  slip  plane  in  the  direc- 
tion x  is  simply  a  particular  case  of  the  expression  of  u  (x  ,x  ,x  ), 
as  obtained  before  for  the  screw  dislocation: 


Au2(x2)  =  u2(x2,-)  -  u2(x2,  --)  =   2u2(X2,_)  .  (149) 


A  suitable  form  for  .Au  (x  )  can  be  obtained  either  from 
Equation  (131)  or  the  set  of  Equations  (133)  to  (136): 


3"4V  an     2 


(x  -qa)   +  a 


2 
(x2-qa) 


(150) 
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Au2(x2)    =    Au°(x2)    +    AuJ(x2)    +    Au^(x2)    +    Au^)    +    .. 


(151) 


with 


.    0  a     I  -1   X2  -1   X2    a^l  a"      1  X2 

Au„(x  )   =  o~  S17  +  tan       —  +  tan       )   +  s~  TT,  ~^ — o 

2      2  2rr   I  a  a       (         8tt  1-V     2      2 


2      2  ,2      2 

1      3-4V    J        .     V3  ,  .    .      (Va)    +S 

1ft?  1=7     X2  *»  —2-  +  (X2"a)  * "IT- 

x2  (x2-a) 


(152) 


,    1,      .  a      3-4V 

Au„(x„) 


2      2  32tt    1-V 


2      2  ,  x2  2' 

X2+a  ^X2~3        +    3 

^7!  — ^—  +  Bn g — 

X2  (X2  "  3) 


a3       1 


16rr  1-v 


2  2    '  ,2  2 

Ux0  +   a  (x  -a)      +   a  J 


(153) 


.    2.  a       3-4V 

AW    =    "96^1=7 


xQ        x0-a  2      2  2      2 

_2  2  x2+a  (xo~a)    +a_ 


a4       1 


48tt  1-V 


l(£.»)s     ^V]  J 


3  a         3-4V 

AVV    =    288^  117 


2        2 
x2(x2-3a  ) 


3    '  3         .2        2    3 

2         (X2"a)  (X2  +  a    } 


(X2-S)    C(X2-a)2"3a2]'\  a6      J_    /VVa2) 

r,     .2  2f      I    +240TT1-V  U  2  2,\4 

l_(x2-a)    +a    J  J  M^x2+a   ^ 

(X2_a)    [(x2-a)2-a2]^ 
[(x2-a)2  +  a2]4       J       ' 


(154) 


(155) 
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The  curves  representing  Au   are  plotted  in  Figures  22  and  23, 
corresponding  to  the  region  close  to  the  dislocation  line  and  to  a  more 

extended  region,  respectively.   Only  a  few  terms  Au   are  needed  to 

_3 

obtain  relative  precision  equal  to  10    in  the  region  outside  the  points 

x„  =  -a  and  x  =  2a.   Unfortunately,  the  region  in  between  is  not  known 
and  Au  (0)  and  Au  (a)  must  be  evaluated  by  an  interpolation  scheme. 
Values  of  Au   are  listed  in  Table  7. 

TABLE  7.   Relative  Displacement  Across  the  Slip  Plane 
for  an  Edge  Dislocation 

X2/a  AU2/a 

-  5  .0512 

-  4  .0626 

-  3  .0805 

-  2  . 1123 

-  1  . 1817 

0  «   .4 

1  «   .68 

2  . 8280 

3  .8921 

4  .9221 

5  .9391 


The  width  of  the  dislocation  can  be  reached  by  evaluating  the 
region  where  the  relative  displacement  has  values  between  a/4  and  3a/4. 
From  Figure  22,  this  region  can  be  easily  measured  and  has  the  value 
W  =  1.93a.   This  value  is  slightly  higher  than  that  obtained  from 
Peierls'  model  (W  =  1.5a  for  V  =  1/3).   This  would  mean  that  our  model 
shows  a  dislocation  slightly  more  extended  than  Peierls'  model.   This 
result  is  opposite  to  what  has  been  found  for  the  screw  dislocation. 
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The  relative  displacement  at  atomic  points  can  also  be  obtained 
directly  from  Equation  (128): 

4.  CO  f    f     CO  CO 

2       +„        j 

q=0         q= 


Au    (na)    =    pa  2  I     +      E      f~G      (pa,na-qa,0) 

p-  _  cc    V  \q=0        q=l/ 

-  G22(pa,na-qa,a)J    -  (_G23 (pa, na  , 0)    +  G23(pa,na,  a)J  I      .  (156) 

Replacing  the  variable  (n-q)  by  u  and  making  the  singular  point  appear 
at  q = n,  leads  to  the  three  following  equations  for  Au  (na)  correspond- 
ing to  n>l,  n=l  and  n  =  0,  respectively. 

2       +m       f  n  n-1  °°  \ 

Au2(na)   =   |ia  E  E+E+E       QG22(pa,ua,0) 

p=  -  °°  Vu=l        u=l        u=l/ 

+  co 

"1  2 

-  G22(pa,ua,a)J    -pa  E        G23(pa,na,a) 

p—    _  CO 

00  2 

+    4pa        E     (^G22(pa,0,0)  -G22(pa,0,a)]   +^--2p^2   G22(0,0,a), 
p=l  2 

(157) 


Au2(a)   =   2  pa  E  E     Q}      (pa,ua,0)    -  G?2  (pa,ua ,  a)j 

p=  -  co  u=i 


2  r  "i 

+    pa  E       |j322(pa,a,0)    -  G22(pa,a,a)    -  G23(pa,a,a)J 


2 
P 


+    4pa2      E     QG22(pa,0,0)    -   G22(pa,0,a7|     +  ^~  -  2pa2   G^CO.O.a), 
p=l  2 

(158) 

_|_CO  CO 

Au2(0)    =   2pa  E  E     Q}      (pa.ua ,0)    -  G22(pa ,ua, a)] 


—     —    CO 


P= 

ro  2 

+   2pa        E     [c22(pa,0,0)  -  G22(pa,0f  a)]  +  J£-  -  pa2   G22(0,0,a) 
p=l  2 


(159) 
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After   computation   of   the  discrete   sums   using  Equation    (26) , 
these  three   equations   become,    in   the   same  order: 


,      ,         2ua  a     11.4392    -   9.4720V         a       1         n 

Au2(na)    =  — 


8tt 


l-v 


8TT   1-V      2    , 
n  +1 


16rr 


E     +        E 


n       p=n+l/ 


l-v     '      2 


l-v     2    , 
P  +1 


(160) 


,    „         2ua  a       6.4060    -   5.6512v 

Au2(a)  =  -ft  +  st  r^ 


(161) 


.      ,„.       ua             a      11.4392   -   9.4720V 
Au2(0)    =—  +  —  — 


(162) 


For  each  value  of  n,  a  direct  comparison  can  be  made  between 
Equations  (160)  to  (162)  on  the  one  hand,  and  Equations  (133)  to  (136) 
on  the  other  hand.   Since  both  ought  to  be  identical,  C   can  be  easily 
deduced  by  subtraction.   Values  of  Au  (na)  computed  in  both  ways  and 
the  corresponding  C  (n)  are  listed  in  Table  8.   These  values  have  been 
computed  for  v  =  1/3.   The  force  constants  C   seem  to  be  much  more  sen- 
sitive to  n  in  the  case  of  the  edge  dislocation  than  in  the  case  of  the 
screw  dislocation   Therefore  the  evaluation  of  the  atomic  displacements 
at  x  =  0  and  x  =  a,  using  the  technique  employed  for  the  screw  dislo- 
cation is  not  possible  due  to  the  uncertainty  of  C_(0)  and  C  (1). 
If  the  variation  of  C  (n)  is  neglected  after  a  few  atomic  distances 
from  the  dislocation  line,  and  C  (n)  set  equal  to  C „(co)  for  every  n, 
the  relative  displacement  becomes 


Au  (na)  =  a  -  -^—  I  E  +   E 
2  16tt  1  p=n   p=n+l. 


3-4v  „  p  +1    1 
m  — =—  + 


l-v 


1-V   2  ^ 

p  +1  J 


a   1 


8rr  1-V   2 

n  +1 


(163) 


Table  8.   Relative  Displacements  and  Force  Constants  C   at 
Singular  Points  for  an  Edge  Dislocation 
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0 
1 
2 
3 
4 
5 
6 
7 
8 
9 
10 


l^oo-Sg. 


253 

232 

3625 

4027 

4243 

4376 

4466 

4531 

4580 

4619 

4649 


Au    (n)/a 

C2 

.4 

sa     6.80 

;       .68 

as     4.40 

.8280 

4.2792 

.8921 

4.0867 

.9221 

4.0176 

.9391 

3.9872 

.9500 

3.9729 

.9576 

3.9635 

.9632 

3.9588 

.9675 

3.9557 

.9709 

3.9526 

4927 


1.0000 


3.9432 


The  relative  displacement  at  atomic  points  corresponding  to 
negative  x   (x  =  -na)  can  be  deduced  directly  from  Equation  (150) : 


£u0(-na) 


16-:  |  S  +    E  J  |  1-, 
(p=n   p=n+l 


3-4V  .  p  +1    1 
m  — s—  + 


1-V   2 

P  +U 


a   1 


8tt  1-V  2 

n  +1 


(164) 


One  can  notice  that  adding  Equations  (163)  and  (164)  leads  simply 
to  a.   This  means  that  making  the  approximation  that  C  (n)  is  constant 
and  equal  to  C  (ro)  for  every  n  would  imply  that  the  relative  displacement 
is  symmetrical  with  respect  to  x  =  0.   But  this  is  in  contradiction  with 
the  actual  relative  displacement  computed  from  Equations  (152)  to  (155) 
and  listed  in  Table  8,  where  a  symmetry  with  respect  to  x  =  a/2  is 
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clearly  apparent.   So,  contrary  to  the  case  of  the  screw  dislocation, 
C„(n)  cannot  be  replaced  by  C  (ro)  ,  and  their  difference  is  sufficient 
to  shift  the  symmetry  of  the  dislocation  with  respect  to  x  =  a/2  to 
x  =  0.   In  fact,  it  seems  logical  that  in  both  cases,  for  the  screw  disloca- 
tion as  well  as  for  the  edge  dislocation,  the  dislocation  line  lies 
between  two  rows  of  atoms. 

As  mentioned  for  the  screw  dislocation,  a  distribution  function 
for  infinitesimal  dislocation  loops  can  be  introduced,  following 
Equation  (64).   Its  expression  is  the  component  o,„1  of  the  dislocation 
density  tensor  mentioned  previously.   So,  by  differentiating  Equations 
(152)  to  (155),  this  distribution  function  a      (x  )  has  the  form 


a21(X2)  =  a21U2)  +  Q'21(X2)  +  Q21(X2}  +  a21(X2)  + 


(165) 


with 


0         1 

Q'21CX2;  -    8tt(1-v) 


2   2     2,      2 
La  -t-x    a  +(x  -a)  _ 


l   3-4V 
16tt  l-V 


Qm 


2   2 

x2+a 


+  im 


2   2. 
(x  -a)  +a 

2~ 
(x2-a) 


o      2   2 
2     a  -x„ 

i    l      2 


8TT  l-V    2      2 

(x2+a  ) 


(166) 


1  ,   .     a   3-4v 
°'21(X2)  =  18n  T37 


2   + 


2   2,     ,2   2 

ux  +a    (x  -a)  +a 


X2~tj 


a 1_ 

8rr  l-V 


2  2   2 

,  2   2^     (x  -a)  +a 
_(x2+a  )       2 


(167) 


2         a   3-4V 
*21(V  _  96tt  T=v" 


2   2 
a  -x„ 


2  ,     ,2    ,  2   2  2 
x2  (x2-a)     (a  +x2) 


2   ,     ^ 
a  -  (x9-a) 

[_a  +(x2-a)  J 
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48rr  1-v 


2   2  3 


(x0-a)' 


(x„+a  )     r  ,2   2"] 

2         L(x2_a)  +a  J  - 


(168) 


a   3-4v 


21 


960tt  1-v 


2    (X2"a) 


„  2  2   4   4 
3a  x  -x  -a 
2   2 

,  2   2  4 
(x2+a  ) 


2       4         4   4' 
3a  (x  -a)   -  (x -a)  -a 


Qx2-a)  +a  J 


2       2         4   4' 
10a  (x  -a)  -5(x  -a)  -a 


Q(x2-a)  +a  J 


r   2  2    4   4 
10x2a  -5x2-a 


240-  1-v 


,  2    2  3 
(x2  +  a  ) 


(169) 


Successive  approximations  of  a?  (x  )  are  plotted  in  Figure  14. 

A  direct  comparison  with  Peierls'  model  cannot  be  achieved 
successfully  because  of  the  lack  of  symmetry  of  the  displacement  field 
around  the  core  in  our  model.   Both  models  give  the  same  result  for 
points  far  from  the  dislocation  line,  but  cannot  be  matched  close  to 
the  core  region.   We  shall  see  that  a  description  of  the   edge  disloca- 
tion from  an  array  of  prismatic  loops  is  much  more  satisfactory  and 
more  close  to  the  real  atomic  arrangement  at  the  center  of  the  defect. 


Self-Energy  of  the  Edge  Dislocation 

As  seen  previously,  two  steps  are  required  to  obtain  the  final 
energy  of  the  edge  dislocation.   First  the  energy  of  the  system  of 
forces  will  be  computed,  and  then  a  correction  energy  term  will  be 
introduced  to  subtract  the  excessive  strain  energy  across  the  slip 
plane,  due  to  the  system  of  forces  itself. 
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When  the  rectangular  array  of  forces  is  extended  in  the  direc- 
tion perpendicular  to  the  Burgers  vector,  it  simulates  a  system  composed 
of  two  antiparallel  straight  edge  dislocations  separated  by  La.   The 
energy  per  unit  length  of  the  array  forces  is  the  limit  of  W/Ra  when 
R  is  much  larger  than  L.   From  Equation  (28) ,  it  takes  the  form, 


W  2 

=  Lpa 


Ra 


Ua 

LC2 


4tt 


|ia 


2rr(l-v) 


£n  L  +    \i& 


2(ia 


7.6170  -  2.8588V 
4TT(1  -  V) 


with 


A  = 


2.860  -  2.368V 
1  -  V 


C2    C3 


(170) 


(171) 


As  previously  seen  for  the  energy  of  the  kink  in  a  screw  dislocation, 
C_  is  an  unknown  parameter  which  cannot  be  obtained  by  a  physical 
argument  as  in  the  case  of  C  . 

The  choice  of  a  suitable  displacement  field  describing  the 
dislocation  will  be  made  as  for  the  evaluation  of  the  correction  energy 
of  the  double  kink  in  a  screw  dislocation.   The  region  where  the  relative 
displacement  across  the  slip  plane  is  larger  than  a/2  is  contained 
between  the  planes  x  =  a/2  and  x  =  La  -  a/2,  as  shown  in  Figure  25. 

A  displacement  field  chosen  to  approximate  the  actual  relative 
displacements  of  atoms  across  the  slip  plane  is: 

(172) 


v1(xi,y2,x3)  =  0, 


wvv  =  w  w  -  v 


v„(x  ,y„,x„)  =  u„(x  ,x„.x„)  , 
3   1,J2'  3      3   1'  2'  3   ' 


(173) 
(174) 
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depending  on  whether  x   is  positive  or  negative,  y   is  taken  as: 


y2  =  X2  +  2 
a 

yo  =  xo  -  o 


for   x  >  0 


for   x  <  0 


(175) 
(176) 


The  principal  difficulty  with  this  choice  of  v  is  that  it  is  not 
symmetric  with  respect  to  x  =  0 ,  since  u   varies  with  x  .   We  have: 


and 


where 


VwV  =  WW 


v2(xl'y2'"X3)    =    -U2(Xl'X2+a'X3)    +   V 


|u2(x1,x2,x3)|  4    |u2(x1>x2+a,x3)|  . 


(177) 


(178) 


(179) 


Since  this  lack  of  symmetry  will  have  an  effect  only  on  the  constant 
terms  which  do  not  appear  in  the  coefficients  of  R  or  £,n   R,  it  will  be 
considered  a  sufficient  approximation  for  our  study. 

The  tensors  a   and  T  corresponding  to  the  displacement  fields 
u  and  v,  respectively,  have  the  following  components: 

f  ,  /S    aU3 


22 


'33 


9u2      Su3 
2        3 


5u      9u 
(9+2^)  ^±  +    X  tt-=- 


>2+    BU3 

dX3+  ^ 


(180) 
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fTll(y3'V  =  CT11(X2'X3)  ' 


T22(y2'X3)  =  CT22(X2'X3)  ' 


33  y2'    3      33   2'  3   ' 


V>T23(y2'X3)  =  a23(X2'X3)  "  * 


(181) 


The  correction  energy  is  by  definition  the  difference  between  the  strain 
energies  associated  with  the  above  stress  fields  contained  within  the 
volume  of  integration,  i.e., 


E  =  i  f  a.  .    e. .  dV  -  ~  f  T. .  e. .  dV1  . 

c   2j   ij   ij      2  ^,   ij   ij 


(182) 


By  using  Green's  theorem,  the  volume  integrals  can  be  transformed  into 
surface  integrals  on  the  surfaces  bounding  V  and  V1 : 


=  i   f  a.  .   u.  n.  dS  -  i  f  T. .  v.  n.  dS* , 
2  «J   ij   i   J      2  J  ,   ij   i   j 


(183) 


where  n  is  the  normal  to  the  surface  of  integration.   The  differences 
between  S  and  S'  are  mainly  due  to  a  translation  of  the  plane  x  =  a/2 
of  an  amount  +a/2,  and  a  translation  of  the  plane  x  =  -a/2  of  an 
amount  -a/2.   Because  of  these  differences  between  the  limits  of 
integration,  all  integrals  over  a  function  y   can  be  transformed  into 
integrals  over  a  function  of  x  : 


La  La-a/2 

x3>0,   /    f(y2)  dy2  =  J       f(x2)  dx2, 
a  a/2 


(184) 


La-a  La-a/2 

X3  <  °'   /    f(y2)  dy2  =  / 

J  0       *  Z   a/2 


f(x2)  dx2 


(185) 


So  the  correction  energy  has  the  form 
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E  =  -  f  (cr  .  u   -  T  .  v. )  n.  dS 
c   2  -g   ij   1    ij   i   J 


(186) 


Replacing  v  and  T  in  function  of  u  and  a,    respectively,  in  Equation 
(186)  leads  to 


E      a/2 
c      C 


Ra 


-a/2 


CX3a22(l'X3)  +  ^Vl'VJ  dXc 


La-a/2 

+  I  CI  CT23(X2'f)  +  ^2(X2'I)  "  *T3dx2> 

a/2 


(187) 


E        a/2 
-£  =  -  2  f 

Ra      o   L 


ou 


(X+2u)  x 


(£,x.)  +  (^i-\)  u„  (|,x„)|  dx. 


3  ^   '2'A3 


3  v2'  3 


La/2 
.      ,a  a       p 

"  X  a  u3(2'2}  +  ^  J 


5U2 

ox    2  2 


U2(X2'2) 


-  1 


dx„ 


+  ^a  (~u (La/2,  a/2)  -  u3(a/2,a/2)J  .  (188) 

For  computing  these  integrals,  we  shall  use  the  expressions  for 
u.  and  its  derivatives  from  Equations  (131)  and  (132) .  The  final  summa- 
tion with  respect  to  q  will  be  completed  at  the  last  step.  Expressions 
used  for  x  (cu  /ox  )  and   Bu  /ox   are  listed  below. 

O        w        i-  ._,        *J 


du2  ax3  /  L 

<3  (ST  (W  =  16*  E n      1 

2  Vq=0    q=l 


x0  -  qa 


^  .    <S2 
(x3-qa)  +(x3+-) 


x„  -qa 


^2  ,     aN2 
(x2-qa)  +(x3--)  ^ 


2 


(x2-qa) 


a,  2 


[(x2-qa)  +(x3-|) 


(x2  -qa) 


_  Qx2-qa)2+(x3+|)2] 


a.  2 


16tt(1-v)  \(X3  ~2} 


l*3~S> 


r   2  ,    <\2T 
LV(X3"2>  J 
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+    (X3+2} 


,        a, 2  2 

(X3+2}       "   X2 


C&'vW 


(189) 


3112    t  1  a 


/  n2       ,  aN.2 

(x2-qa)     +(x3--) 


2 
1-V 


/•  n2    /  S2 

(x2-qa)    +(x3  +  -) 


.  av.3 

(X3+2) 


[(x   -qa)2+(x  +|)2J 


r  aA2 

(VJJ 


16tt(1-V) 


[(x2-qa)2+(x3-|)23 

r  2  3x2-.  r  2        /  a.  2 -j 

•2   LX2-(x3-5)     J         X2LX2    "    (X3  +  2)     J 


r  2  a  2~]  f  2  a,2n: 


(190) 


We  will   not   reproduce  here   the   details   of    integration,    but   only  write 
the   final    results   for   each   step   of   the   computation.      The   correction 
energy  becomes 

'  L 


E  2 

c  _    ua 

Ra   ~     4tt 


(q-i)2+l 


n   +         J^   (q"2)^- O 

Q=0         q=l/  (q-o> 


+   4   tan"      (q  --)} 


ua      1-2V 
2tt     1-V 


Sm  L 


p« 


Ua 


2      L  +     2n 


1       Z1  «    5        7 
1^     2  ^  4  +   40 


,      -1    _        1    16v     -    18V  +5  V       „    „ 

+   tt  -   tan        2   +  —  — -z-rr. „.  ,      +  ■  #72  2 


4  ( 1  -  v)  ( 1  -  2  v) 


1-2V 


J 


(191) 


Computing  the  single  sums  as  mentioned  in  Equation  (26)  leads  to  the 
final  expression  for  the  correction  energy  of  the  system,  per  unit 
length  of  edge  dislocation. 


90 


Ec         ua2L         i^a2      1      .  |ja2     r  .  5198   +    2.3167V  -    .3182V' 

^ L  +  -w  L— 


Ra 


rr     l-v~"~   '     2rr     L  (l-v)(l-2v) 

Recalling   Equation    (170),    the   total    energy  per   unit   length  becomes 

Ra        Ra        Ra    ' 


J  .         (192) 


(193) 


that   is, 

E 


Ra 


T  2 

=   L   \xa 


\i&         A         1 
C2   +   4TT  "   2 


2      , 
ua        1 

2n    1-v 


&n  L 


2rr 


2  r 


4tt    pa    |i  -  JL 
2  3 


4.3273    -   6.7297V  +    2.5406V 
(1-V)(1-2V) 


2-1 


(194) 


For  E  /Ra  to  represent  the  energy  of  the  system,  C  has  to  be 
T  « 

chosen  such  that  the  term  linear  in  L  vanishes.   Under  this  condition, 
C     takes  the  same  value  as  for  the  screw  dislocation.   Unfortunately, 
C„  remains  unknown.   It  is  clear  here  that  it  cannot  be  set  equal  to 
C  as  it  is  a  priori  assumed  for  the  energy  of  the  kink  in  the  screw 
dislocation.   If  it  were,  the  core  parameter  would  be  unrealistically 
large  if  one  expects  the  same  order  of  magnitude  as  for  the  screw  dis- 
location.  We  can  even  foresee  that  C  has  to  take  a  value  much  larger 
than  C  ;  this  force  constant  C  ,  computed  from  the  displacement  field, 
was  increasing  greatly  for  atoms  close  to  the  dislocation  line.   Since 
C   corresponds  to  the  closest  atoms  to  the  dislocation  line,  it  can  be 
expected  to  be  of  the  order  of  magnitude  of  C  (0) . 

However,  the  total  energy  of  the  system  composed  of  two  parallel 
edge  dislocations,  and  of  opposite  signs,  has  the  form  expected. 


(ia 


7   La 


Ra   2n(l-V)    r 


(195) 
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and  the  self-energy  for  a  pure  edge  dislocation  is  simply 

!g  =    ^    r(k  —   .  (196) 

Ra   4tt(1-v)    rQ 

The  fact  that  the  essential  features  of  the  self-energy  of  the  edge 
dislocation  are  reached  with  our  model  is  very  encouraging.   The 
evaluation  of  the  force  constant  C  remains  the  only  uncertainty  in  the 
energy  of  the  edge  dislocation. 

Single  and  Double  Kinks  in 
an  Edge  Dislocation 

Similar  to  kinks  in  a  screw  dislocation,  a  single  kink  or  a 

double  kink  in  an  edge  dislocation  can  be  simulated  by  adding  only  one 

shear  loop  at  the  end,  x  =0,  of  each  row  in  the  region  x  <  0  for 

1  ^ 

a  single  kink  and  in  the  region  -Na  <  x  <  Na  for  a  double  kink. 
The  array  of  forces  takes  the  configuration  shown  in  Figures  26  and  27, 
corresponding  to  single  and  double  kinks,  respectively. 
A.   Displacement  field 

The  displacement  field  corresponding  to  these  defects  is  given 
by  the  superposition  of  the  displacement  field  of  the  edge  dislocation 
and  the  displacement  field  of  the  extra  forces  introduced  to  simulate 
the  defect.   The  latter,  u' ,  has  the  following  expression,  for  a  single 
kink  or  a  double  kink  of  length  2Na,  respectively. 

U1(SK)  =  ^  Zn  f(p)  (197) 

p=0 

and 

u!(DK)  =  ±f-  I      f(p)  ,  (198) 

p=-N 
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^fe^S) 


Figure  26.   Array  of  Forces  for  Single  Kink  in  an  Edge  Dislocation 


-^j) — -£y-£3 


~<T~(F 


Figure  27.   Array  of  Forces  for  Double  Kink  in  an  Edge  Dislocation 


where 


f(p)  =  G^x^pa.x^a.Xg--)  -  G.^-pa^+a^  +  -) 


+  G„„(x  -pa,x„ ,x„  -  — )  -  G„„(x  -pa,x„,x„+— ) 
22   1  F  '  2'  3   2      22   1  K  '  2'  3   2 


G23(x1-palx2+a,x3-2)  -  G^C^-pa, Va,x3+-) 
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+  G23(x1-paIx2>x3-j)  +  G23(VPa'X2'X3  +  2) 


(199) 


We  shall  only  consider  the  displacements  in  the  x   direction, 
since  we  are  mainly  interested  in  the  change  of  the  atomic  configura- 
tion in  a  plane  parallel  to  the  slip  plane.   Due  to  the  symmetry  of 
the  extra  rows  of  forces  with  respect  to  x  =  -a/2,  we  shall  only  take 
account  of  points  corresponding  to  x  >  -a/2,  for  example.   The  expres- 
sions for  the  displacement  of  points  of  application  of  point  forces 
will  be  different  from  those  for  other  atomic  point,  since  the  former 
are  singular  points.   The  following  equations  obtain  in  the  plane 
x  =  a/2  for  single  and  double  kinks,  respectively. 

(a)   All  x  except  x  =  0  and  x  <  0 


U2(SK)    "   32tt 

u=n 


3-4V 


1-v 


/2~2       ;  ~2~  /2~2~  ~2       2 

a/u   a     +(x  +a)  a/u   a    (x0+a)    +a 


/"2~2      2  /2~2      2      2 

a/u  a  +x  a/u  a  +x  +a 


1 

1-V 


(x2  +    a) 


f  2    2    i  ^1 

|u    a   +(x  +a)     J 


3/2 


(x2+a) (x2+2a) 


x2(x2-a) 


3/2 


3/2 


r  2   2     2  1    '  r   2    2    ,  v2      2T'  r  2   2      2      2  ~| 

lu   a  +x    J  u   a   +(x  +a)    +a    J  [u   a   +x  +a    J 


3/2 
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(b)      x    =  0      and     x     <  0 


V.  Vu  +1  Vu  +2 


l-v 


(u2+l) 


3/2 


-      I    3"4v     I  1 


2 

2C„    ' 


(201) 


and   for    a  double  kink  of   length   2Na, 


(a)      All   x     except   x    =  0      and   -Na  <  x     <  Na 


a  a2       g      I    3-4v 

U2(DK)CX1'X2'2;         32rr  „S    (l-v) 

p=-N    I 


/  2  2 

/(x1~pa)"+(x2+a) 


y'  2  2      2  /  2      2  / 

(x  -pa)    +(x  +a)    +a  (/(x  -pa)    +x         V(x  -pa) 


2      2      2 


(l-v) 


(x2+a) 


3/2 


r  2  2*1  f  2221 

_[_(x   -pa)    +(x2+a)     J  ^(x^^-pa)    +x2)     J 


,3/2 


(x2+a) (x2+2a) 


x2(x2-a) 


r  9   ,  2     2-i3/2         f.  .2     2     2-|3/2 

|^(x   -pa)    +(x  +a)    +a    J  |_(x   -pa)    +x2+a    J        J 


(202) 


(b)      x    =  0      and      -Na  <  x     <  Na 


«    aN  a  2        J        3"4v  ! 

U2(nc)(na'°»2)  a55S     \T\~~T= 

p=-N      i 


1 

l-v 


"1 


V(n-p)    +1 
N+n        N-n\ 


Q(n-p)2+l] 


3/2 


r       2     2-i3/2 
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The  force  constant  C  used  here  ought  to  be  C  (0) ,  computed 
from  Equation  (162).   It  can  only  be  approximated,  since  Au  (0)  for  the 
edge  dislocation  is  only  known  approximately  as  0.4a.   In  any  event,  we 
shall  use  this  value  to  obtain  the  general  trend  of  the  atomic  config- 
uration in  the  immediate  vicinity  of  the  kink.   The  atomic  displacements 
are  tabulated  in  Table  9  for  the  single  kink,  and  in  Tables  10  to  12 
for  the  double  kinks. 


TABLE  9.   Atomic  Displacements  for  a  Single  Kink 
in  an  Edge  Dislocation 


xl 

U2(SK)(W 

x2,0 

x2  =  -a 

x2=-2a 

-  5 

.  1243 

.0854 

.471 

-  4 

.  1242 

.0852 

.0468 

-  3 

.1239 

.0850 

.0464 

-  2 

.1224 

.0834 

.0442 

-  1 

.  1145 

.0764 

.0355 

0 

.0622 

.0428 

.0223 

1 

.0090 

.0090 

.0116 

2 

.0020 

.0020 

.0029 

3 

.0004 

.0004 

.0007 

4 

.0002 

.0002 

.0003 

5 

.0001 

.0001 

.0001 
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TABLE  10.   Atomic  Displacements  for  a  Double  Kink  of  Length  2a 
in  an  Edge  Dislocation 


U2 

(E*0(X1 

V 

x  /a 

x2  -  -2a 

X2  «  ~ 

a 

X2  =  ° 

x2  -  a 

0 

.0404 

.0804 

.  1194 

.0404 

1 

.0338 

.0719 

.1109 

.0338 

2 

.0115 

.0120 

.0120 

.0115 

3 

.0031 

.0022 

.0022 

.0031 

4 

.0010 

.0006 

.0006 

.0010 

5 

.0003 

.0002 

.0002 

.0003 

TABLE  11.   Atomic  Displacements  for  a  Double  Kink  of  Length  4a 
in  an  Edge  Dislocation 


U2(DK)(X1 

x2)/a 

x  /a 

x2  =  -2a 

x2  =  -a 

X2  =  ° 

x2  -  a 

0 

.0446 

.0835. 

.  1225 

.0446 

1 

.  0432 

.0824 

.1214 

.0432 

2 

.0348 

.0725 

.  1115 

.0348 

3 

.0118 

.0122 

.0122 

.0118 

4 

.0031 

.0022 

.0022 

.0031 

5 

.0010 

.0006 

.0006 

.0010 

TABLE  12.   Atomic  Displacements  for  a  Double  Kink  of  Length  6a 
in  an  Edge  Dislocation 


U2(DK)(X1 

x2)/a 

x  /a 

x2  =  -2a 

x2  =  -a 

X2  =  ° 

x2  =  a 

0 

.0460 

.0844 

.  1234 

.0460 

1 

.0456 

.0842 

.  1232 

.0456 

2 

.0435 

.0827 

.  1217 

.0435 

3 

.0348 

.0743 

.1133 

.0348 

4 

.0118 

.0123 

.0123 

.0118 

5 

.0031 

.0023 

.0023 

.0031 

97 


Figures  28  to  31  show  the  atomic  configurations  of  the  defects 
in  the  planes  x  =  ±  a/2,  corresponding  to  the  values  listed  in 
Tables  9  to  12  in  the  same  order.   The  same  comments  made  for  kinks 
in  the  screw  dislocation  apply  here,  that  is,  particularly  the  local- 
ization of  the  distortion  in  the  immediate  vicinity  of  the  defect. 

B.   Energy  of  the  double  and  single  kinks 

The  energy  of  a  double  kink  will  be  defined  as  the  difference 
between  the  energy  of  the  modified  rectangular  array  of  forces  and  the 
unmodified  rectangular  dislocation  loop.   First,  the  energy  difference 
between  the  two  arrays  of  forces  is  equal  to  the  self -energy  of  the 
extra  double  row  of  forces  added  to  their  interaction  energies  with  the 
rest  of  the  array.   Their  expressions  are 


2  4 
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(206) 


In   the   expressions   above,    we   always   considered   the  length   of 
the  double  kink   to  be  much   smaller   than   the  length   of   the  dislocation. 
Replacing   the  Green's   functions   by   their   explicit    expressions   leads   to 
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(P   +D         / 


(207) 


The   term    pa/C    ,    depending  on  the   force   constant,    can  be 
replaced  by   its   value,    as   defined   in  Equations    (113)    and    (114). 
The  final    expression  for   the   energy   difference  between  the   systems 
of   forces  becomes 

3  I  R/2 
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1-v 


Jv  +< 


103 


+  1-V  I  ,  2    3/2     2   3/2 
(P  +2)       (p  +1) 


Ita3  2g+1    3-4V  I        p      ,\  +   1    (    2P 

(208) 

If  N  becomes  large,  but  still  smaller  than  the  length   of  the  edge 
dislocation,  the  double  kink  behaves  like  two  separate  single  kinks. 
The  corresponding  energy  of  the  system  of  forces  for  each  single  kink 
is  the  limit  of  W/2  when  N  increases: 


+  T^  l7T^37g-   2  >2'^ 

\(P  +2)       (p  +1) 

The  interaction  energy  between  the  two  kinks  is  simply  (2W   -  W) . 

As  for  the  kinks  in  a  screw  dislocation,  a  correction  energy 
term  has  to  be  introduced  to  remove  the  excess  strain  energy  contained 
in  the  region  between  the  planes  of  forces.   Unfortunately,  the  same 
difficulties  arise  as  before  and  we  shall  not  be  able  to  perform  a 
complete  analysis  of  the  energy  of  a  double  kink.   The  computation  of 
the  strain  energy  can  be  done  exactly  in  the  same  way  as  for  the  screw 
dislocation,  and  its  development  here  would  only  be  a  repetition  of  the 
treatment  of  the  correction  energy  for  a  double  kink  in  a  screw 
dislocation. 

We  anticipate  that  in  the  complete  expression  for  the  energy  of 
the  kink,  the  coefficient  of  N  can  be  made  to  vanish,  leaving  only  a 
pure  number.   It  is  interesting  to  point  out  that  Equations  (208)  and 
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(209)  have  exactly  the  same  form  as  Equations  (115)  and  (116) ,  except 
for  the  presence  of  the  force  constant  C  .   This  is  due  mainly  to  the 
screw  dislocation  character  of  the  kink  segments.   If  a  complete  treat- 
ment of  the  correction  energy  can  be  achieved,  the  energy  of  a  kink  in 
an  edge  dislocation  is  thoroughly  determined. 

The  results  already  obtained  are  very  encouraging  and  a  complete 
treatment  should  lead  to  a  better  understanding  of  kinks  in  edge 
dislocations. 


CHAPTER  6 


EDGE  DISLOCATION  IN  SIMPLE  CUBIC  CRYSTAL  CONSTRUCTED 
FROM  AN  ARRAY  OF  PRISMATIC  LOOPS 


A  rectangular  dislocation  loop  formed  by  four  segments  of  edge 
dislocations  can  be  simulated  by  an  array  of  primitive  prismatic  loops 
as  defined  in  the  first  chapter.   This  is  equivalent  to  a  rectangular 
plane  of  vacancies  plus  an  extra  contraction  in  a  direction  normal  to 
this  plane.   The  loop  is  shown  in  Figures  32  and  33.   The  forces  F  and 
G  are  those  computed  from  Equations  (15)  and  (16) ,  and  have  the  resepc- 
tive  values: 

F  =  — ^  (210) 


and 


G  =  ^-Z   a   .  (211) 


The  same  procedure  as  used  in  the  previous  chapters  will  be  followed 
here  to  obtain  the  displacement  field,  width  and  self-energy  of  the 
pure  edge  dislocation. 

Displacement  Field 

The  loop  considered  will  be  rectangular,  of  dimensions  La  and 
Ra  in  directions  parallel  to  x  and  x   respectively.   This  means  that 
the  loop  is  composed  of  (R+l) (L+l)  vacancies.   The  edge  dislocation 
considered  will  have  its  line  parallel  to  the  x   axis  and  located 
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Figure  32.   Array  of  Prismatic  Loop  in  x  =  0  Plane 
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Figure  33.   Array  of  Prismatic  Loop  in  x.  =0  Plane 
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near  x  =  0.   So  the  coordinates  x   and  x  will  always  be  considered 
small  with  respect  to  La  and  Ra ,  respectively. 

Considering  the  displacement  field  in  this  region  as  the 
superposition  of  the  displacement  fields  of  all  the  point  forces  leads 
to  the  following  formal  expression  in  terms  of  Green's  functions: 

Xa2      L/2  f  a 

WW  =  -r    E   ,„    "  G3i(xi'x2-qa-x3-2) 

q=  -  L/2    ^     . 

-    G3i(Xl'X2-qa'X3  +  f)    +    G3i(Xl'X2-qa'X3+l+Ra) 

3a  1         Xa2   R+1    r  la  a 

+    G3i(Xl'X2-qa'X3  +  — +Ra)j   +   —  ^  Ki(Xl'X2+i-'X3-2  +  Ra) 

+    G2i(Xl'X2  +  ¥+a'X3-f+ra)    "   G2i(XrX2-¥'X3-f+ra) 

La  a  -\         X+2u     2   R+}        L/y2 

J  r=l      q=  -  L/2 

G  ,(x  +a,x  -qa,x    -~+ra)    -  G     (x  -a,x  -qa,x    --  + ra)  )   ,    (212) 

where  G. . (X  ,X  , X  )  is  defined  by  Equation  (3). 
ij   X  2   3 

Again,  Euler's  formula,  Equation  (29),  will  be  used  to  evaluate 
the  discrete  summations.   The  first  summation  considered  will  be  along 
the  line  of  the  dislocation.   This  sum  is  found  to  be  identical  to  the 
first  term  of  Euler's  approximation: 

L/2  L/2 

Z     f(q)  =  J  f (x)  dx  ,  (213) 

q=  -  L/2        -  L/2 

so  an  expression  of  u.(x  x  , x  )  can  be  written  as: 
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u  (x  ,x  ,x  )  =  -  

1    '   '        8tt(1-v)  (1-2V) 


Xl(X3-2} 


Xl(X3  +  2} 


2  ,  a->2     2  /     ax2 

V'^"^     X1+(X3  +  2) 


R+l 


fa- 


<  n2  /      a     ^2 

,3-4V    (xra)  +(x3--+ra) 

f  16tt   A   l-2v  *"  .     ,2      a    ,2 
'I  (x  +a)  +(x3_2  +  >"a^ 


1-2V 


(x  +a)' 


(xi-a)' 


/•     x2  /     a     n2  n2  ,     a     x2 

(x^a)  +(x3~-+ra)     (x  -a)  +  (x3~-+ra) 


(214) 


U2(X1'X2'X3)  =  °« 


(215) 


u„(x  ,x„,x„)  =  - 

3V  1'  2'  3J        16n(l-V) 


fo 


/~2  .    2 


3-4V 
1-2V 


1-2V 


L/2  +  VL  /4+R 


a. 2 


a.  2 


4^R+^.[x1+(x3-|)  ]  +  Bn    [xi+(x3  +  |)  ] 


2      a,2    2      a.2 

Xl+(X3-2}     V(X3  +  2) 


R+l  f(x3  --+  ra)  (x  +a) 


(x3--+ra)(xra) 


8n(l-2v)     \      a     2       2        a2       2/ 
r=l  i(x3-2  +  ra>  +(x  +a)     (x3~2  +  ra^  +(xi~a)  I 

(216) 

The  second  step  of  the  computation  consists  of  applying  Euler's 
formula  a  second  time,  giving  a  series  of  expressions  which  represent 
successively  higher  approximations  to  the  displacement  field  of  the 
edge  dislocation.   The  mathematical  difficulties  encountered  when  com- 
puting integrals  over  a  region  containing  singular  points  are  solved 
in  the  same  manner  as  for  the  screw  dislocation  in  Chapter  3.   The 


final  analytic  expressions  for  u   and  u   take  the  form 
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U     =   U     +   U     +    U     +    U     +    .  .  .     ,  (217) 


with 


WW  =  is<  h*i +  (va)  tan    x7n  -  (xra)  tan    ^T 

3-4V      ,  a,  (Va)     +(X3+2} 

(X„  +-)    S.T, 


16n(l-2v)         3      2'  ""   ,  x2      ,  ax2 

(xra)     +  (x3  +  -) 


'        Wl*  X1(X3  +  I) 


Jrr(l-V)  (1-2V)     A     2    ,  a,  2  2     .  a, 2 

l^xi+(x3--)  V(X3  +  2^ 


(218) 


The   terms   -to,    -ttx     and  rra    in  Equation    (218)    correspond   to   the  regions 
x     >  a,    -a  <  x     <  a   and  x     <  -a,    respectively. 

1,  y  a(3-4v)  (Va)     +(X3+2} 

WW    =    -   32tt(1-2v^: ^— T2 

(xra)     +  (x3  +  -) 


2  2 

(x  -a)  (x  +a) 


16tt(1-2v)     )  2  a   2  2  a  2[    ' 

^(xra)     +(x3  +  -)  (x1+a)     +(x3+2)"j 


2(  ,  a2(3-4.)  J  Va/2  Va/2 

u, (x1 ,x_,x_)    =   - 


(219) 


lv    1'    2'    3'  96n(l-2v)         \    ,  ,2  a,2  ,2  a   2 

I    (x  -a)    +(x3+2^  (x  +a)    +(x3  +  ^ 

2  f  (x   -a)2 

a  ,  a->  /  1 

(x0  +  -)  <  jr- 


48tt(1-2v)         3      2      j  2  a   2  "I' 

L[(Va)    +(X3  +  2)     J 


(x+a)2  *) 
= rr>       ,  (220) 


[_(x1+a)2+(x3+|)2] 
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and 


3,  ,  a4(3-4v)         ,  a     f    3(-Va)2-(X3  +  f)2 

U      (X      ,X      ,X     )     =     -  — ; — —     (X      +  — ) 


1      1'    2'    3  2880tt(1-2v)         3     2     \  n  n      3 

a^  2 


L(xra)  +(x3  +  l}  3 


3(x  +a)    -(xo  +  2^ 


Q(x1+a)   +(x3  +  |)     ] 

V)2[(Va)2-(vi)2] 


(x0  +  o) 


240tt(1-2v        3      2'      . 

LL<V8    )+(X3  +  2}     J 


2  r  2  a  2-i 

(x1+a)      [_(x1+a)    -(x3  +  -)    J 

Qcx^a)    +(x3  +  |-)    J 


(221) 


The  other  non-zero   component   of    the  displacement   field   is   written 

similarly: 

0  12  3 

U3  =    U3   +    U3   +    U3   +    U3   +    • * •     '  (222) 


with 


0.  ,  1-2V  .  RL 

WW  =  i^T^)aM 


/2  2 


L/2+JL   /4+R 


f        2  a, 2  2  a,2* 

v(3-4v)a  I        V(X3-2}  „    Xl+U3  +  2) 

#72   S +   #7!  ^ 


16tt(1-2v)(1-V) 


2  2 

Cl  Xl 


8tt(1-2v)  (1-v)     1     2    ,  ax2  2    ,  aN2 

LV(X3-2)  Xl+(X3  +  2), 


r  2  a    2 

(x  -a)    +(x    +-) 

(x  -a)  0n  — 1 o 


16tt(1-2v)    V    1 

I  a 

,  *2    ,  S2^ 

(x  +a)    +(x    +^)       I 

-    (xi+a)   ^ 1 V        ,  (223) 


Ill 


1  x,+a 

1,  v  a  .  a.    )  1 

U„(X         X         X     )      =     — - — —      (X       +— ) 


3      1'    2'    3'         16tt(1-2v)         3      2'  \  ,  .2     .  a   2 

(Ya)    +(x3  +  2) 


(xra)    +(Xg  +  -) 


(224) 


2 
WW    -         96tt(1-2v) 


t2  [(x^a)   Q(xi+a)    -(x3  +  f)    ] 


[(xi+a)2+(x3  +  |)2] 

(xra)[(x1-a)2-(x3+|)23">> 
[(xra)2+(x3  +  |)2]2J 


(225) 


and 


^  n    r,  >4   „,  2,  a.  2  ,       a,4l 

3  a4  |(x+a)    Q(x+a)    -6(x+a)     (Xg  +  ?)  +0X3+-)    J 

u    (x    ,  x    , X    )    = 


3      l'    2'    3'  960tt(1-2V)     \  _  _  -.  4 

(x1+a)    +(x3  +  2)    J 

(va)C(va)4-6(va)2(x3+l)2+(x3+l)4J>l 

-   a /       •  (226) 


[(x1-a)2+(x3+|)2] 


The  consistency  of  this  displacement  field  with  results  already- 
obtained  from  the  ordinary  continuum  model  is  checked  by  considering 
x  and  x  to  be  far  from  the  dislocation  line,  and  by  expanding  u   and  u 
in  a  Taylor  series, 

f(Yei-VV  =  f(xrx3)  +  ei  o^  (xrx3)  +  e3  -&q  ^i'V-       (227) 

where  s   and  e   are  small  with  respect  to  x  and  x   respectively. 
The  following  expressions  are  found: 


WW    =    2r7 
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-TT/2  "\ 

-1  X3  1  X1X3 

+    tan  —  >  +   ■     ,*      .    -4~    ,  (228) 

/0  x     (       4rr(l-vT      2      2' 

TT/2  1  j  xi+x3 


where   -tt/2   and  tt/2   correspond   to   x     >  0   and  x     <  0,    respectively. 


Also 


U2(X1'X2'X3)    =   °    *  (229) 


W  W=  stTTT^T  a  {  2  0n 


L/24VW4+R2 


2 


-^(VXI>?     -  4!^  "2^2    ■  (230) 

J  VX3 

These  equations  match  exactly  with  the  ordinary  continuum  model, 
Equations  (146),  (147)  and  (148),  when  one  realizes  that  u  (x  ,x  ,0) 
is  equal  to  -a/4  from  the  way  the  displacement  field  has  been  defined. 

By  contrast  with  the  simulation  of  an  edge  dislocation  by  an 
array  of  shear  loops,  in  the  present  case  the  displacement  field  is  . 
symmetrical  with  respect  to  the  plane  x  =  0.   The  only  singular  points 
occur  for  x  =  ±  a   and  x  =  -na  +  a/2  (n  >  1) ,  and  x  =  0  with  x  =  a/2. 
We  shall  demonstrate  that  the  plane  x  =  0  plays  a  unique  role  for  the 
dislocation  in  the  sense  that  the  displacement  of  singular  points  will 
be  shown  from  the  displacement  of  the  points  symmetric   with  respect  to 
x  =  0.   We  have  the  fundamental  relations: 

u1(x1'x2'na+  2}    +   Ul(Kl'X2~na  ~f^  =  ~2'  (231) 


and 


u3(x1,x2,na  +  -)  =  u3(x1,x2,-na --).  (232) 


Let  us  first  consider  the  special  case  of  n  =  0  for  u  . 
From  Equation  (214) ,  we  have 
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.a.         va         1     a   ^ 

VW^  =  8tt(i-*i-2v)  "2~T  +  T&R    E  f(ra)' 

x  +a        r=l 


(233) 


ax         0° 

VX  'X  ■  -§>  =  8tt(1-v)(1-2v)  T-2  +  ife  Sf(ra-a), 

x+a        r=l 


(234) 


with 


,     ,2   2  2 
*t      ,    3-4Vr  (Xra)  +r  a     2 

f(ra)  =  T^fe" ^-T2  +  T=S 

(x+a)  +r  a 


(X  +  a) 


,  N2   2  2 

(x+a)  +r  a 


(xx  -a) 


2   2  2 
(x  -a)  +r  a 


(235) 


Noticing  that 


E  f(ra-a)  =   E  f(ra)  +  f (0) , 
r=l  r=l 


(236) 


the  two  identities  follow: 


OO  (-K      — ?>1 

,                a                                     a,          a       _    _       x         a      3-4v         l    1      ' 
VVVS*    +  ul(xl,x2,  --)    =_     Ef(ra)+— _<fc- -j 

r=l  (x„+a) 


(237) 


u-,  (x,  -xo>  -o>-  u-,  Cx,  ,x0>-)  = 


a      3-4v„     Ura) 


072 


1      1'    2'       2  1      1'    2'2         16rr  1-2 V*"'  2 

(x^a) 


(238) 


On  another  hand,  from  Equations  (218)  to  (221) ,  u  (x  ,x  , 


2> 


has  the  exact  value 
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-a/4                 ax, 
.  a  _  \         ,.        va 1_ 

VW  ~2}    -     \    Xl/4  +  8tt(1-V)  (1-2V)   2   2 

/A  X1+a 

a/4  1 


(3-4  V)  a  „  (Va) 

+  ^   ,,  ^..,  071   


32tt(1-2\»  ""  ,     ,2  ' 
(x  +a) 


(239) 


where  the  terms  -a/4,  -x  /4  and  a/4  correspond  to  the  region  x  >  a, 

-a  <  x  <  a,  and  x  <  -a,  respectively.   Comparison  with  Equations  (237) 

and  (238)  leads  to 


u  (x  ,  x„ , )  +  u  (x  ,x„,— )  -    -x  /2 

1   1'  2'  2J  V    1'  2'2         1 


-a/2 
a/2 


(240) 


Accidentally,  an  interesting  identity  has  been  proved,  that  is, 

2 

CO 

1 


1    1   3-4v„  (Va) 


2   I6rr  1-2 v  ""'     N2   8n 

(xx~a)       r=l 


Z  f(ra)  . 


(241) 


This  result  can  easily  be  generalized  to  x  =  na  +  a/2  and  x  =-(na  +  — ) 
for  any  value  of  n.   The  displacements  in  the  x  direction  for  these 


points  have  the  form 
u1(xlfx2,na  +  |)  = 


Va 


x    na 


8TT(1-V)(1-2V)    \      2      2   2 

I  x  +n  a 


x    (na+a) 

~~ 2  2 

x  + (na+a) 


+    16*  E      "      EJ       f(ra)     ' 

\r=  1        r=  1  / 


(242) 


when   it    is   noticed   that 


E     f(ra+na)    =      E     f(ra)    -      E     f (ra) , 
r=l  r=l  r=l 


(243) 


and 
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Ul(Xl'X2''na-2) 


8tt(1-v)(1-2v) 


2   2  2 
x  +n  a 


x  (na+a) 
x  (na+a) 


16tT 


E  +   E   f(ra) 
r=l   r=0/ 


(244) 


So,  by  adding  Equations  (242)  and  (244),  the  same  expression 
as  Equation  (237)  is  found,  which  leads  to 


u1(xi,x2,na  +  -)  +  Ui(xi(x2,-na--) 


-a/2 
a/2 


(245) 


The  same  line  of   reasoning    applies   to  u    ,    leading   to 


u3(X;L,x2,na  +  -)    =   u3(xi,x2,-na--) 


(246) 


Equation  (245)  will  help  us  to  find  the  atomic  displacements 
of  the  singular  points  and  the  force  constant  C  .   The  latter  can  be 
obtained  by  expressing  the  displacement  of  these  singular  points  in 
terms  of  Green's  tensor: 

\a2   L/2    , 
u  (a,-na--)= —    E    (  G31(a,x2-qa,  -na-a) 

q=  -  L/2 


R+l   L/2 


+  G   (a,x  -qa,-naA  +   -~  a   E    E  (g      (2a, x  -qa, -na-a+ra) 
31    2        J  2  r=l  -L/2  V11  l 


G   (0,x  -qa, -na-a+ra)  ) 
11   '  2   '        J 


(247) 


After  summing  on  q,  Equation  (247)  becomes 

n+1 


u^a.-na-f)  =  --^-J^—- 


2       2 
.(n+1)  +1   n+1. 
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\+2\i   a 
2   C 


(248) 


By  comparing  the  values  of  u  (a, -na, -a/2)  from  Equations  (217)  and 
(248) ,  C  can  be  computed  for  each  value  of  n.  Several  values  are 
listed  in  Table  13. 

TABLE  13.   Displacements  and  Force  Constants  C^  at  Sin- 
gular Points  for  an  Edge  Dislocation 


u  (-  —  na)/a 


ct(v 


1/3) 


4100 

4511 

4731 

4833 

4882 

4910 

,4928 

4940 

.4948 

,4954 

5000 


4.692 
4.688 
4.679 
4.674 
4.672 
4.671 
4.669 
4.669 
4.669 
4.669 
4.666 


The  relative  displacement   across    the   slip  plane   can  be 
easily   deduced   and  has    the  form 


AU1(X1}   =  W2}    "  W 


) 


^l'-^ 


(2<*9) 
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Recalling  Equation  (239) ,  Au  has  the  final  form 

9 
a  r       >     a  3~^r      (Xl+a)      a        V        aXl 

Aui(xi}  =  TfrFl=2vfc- ^-4^  (i-v)(i-2v)  T-T  •  (250) 

(x  -a)  x  +a 

Values  of  Au   are  plotted  in  Figure  34.   Finally,  a  mapping  on  the 
x  =  0  plane  of  the  atomic  arrangement  is  shown  in  Figure  35. 

A  direct  comparison  with  Peierls'  model  can  be  attempted, 
since  both  models  are  symmetrical  about  x  =  0.   From  Peierls'  model, 
the  relative  displacement  has  the  form 

a    a     -1  X1 
k(xj  =  -  -  -  tan  X  -J-    ,  (251) 

1   1     2    TT         i 

where 

§  =  2T^r  •  (252) 

We  shall  expand  Equations  (249)  and  (250)  in  powers  of  1,'x   and  shall 
select  §  so  that  both  expressions  match  at  large  values  of  x  .   So 
first  order  this  gives 

a  3-2 v  a 

AVV  ^  — r;-  <253> 

and  for   the  Peierls'    model 

a     5 
Au    (x    )    =   -  -2-   .  (254) 

11  IT    X 

By  matching  Equations  (252)  and  (253)  Peierls'  model  is  identical  to 
ours  for  large  x  ,  if 

In  contrast  with  the  screw  dislocation,  this  is  not  the  value  found  by 
Peierls.   Peierls'  expressions  for  §  =  a/2(l-v)  and  §=  (a/4)  (3-2V)  /  (1-V) 
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are  plotted  in  Figure  34.  For  V  =  1/3,  the  values  of  §  are  .75  and 
.875,  respectively.  So  our  model  for  the  edge  dislocation  is  wider 
than  Peierls'  model. 

It  is  impossible  to  say  which  value  is  the  actual  one,  since 
both  models  are  completely  different  and  involve  different  approxima- 
tions.  As  a  conclusion  we  would  say  that  simulating  edge  dislocation 
by  an  array  of  prismatic  loops  gives  a  displacement  field  which 
exhibits  the  expected  symmetry  relative  to  the  extra  half-plane,  in 
contrast  to  the  result  for  the  shear  loops. 

Self-Energy  of  the  Edge  Dislocation 

The  energy  of  the  system  of  point  forces  is  composed  as  usual 
of  the  sum  of  the  self-energies  of  the  point  forces  less  the  sum  of 
their  pairwise  interaction  energies.   Because  of  the  formidable  length 
of  the  expression,  we  shall  divide  the  whole  equation  into  smaller 
subdivisions. 

(a)   Energy  of  the  array  of  G  forces: 
1 


W_  =  (R+1)(L+1)G 
G 


-  G   (2a, 0,0) 
LI 


L 

2 
q=l 

R 

E 

q=l 

R   L 


+  2G  (R+l)  Z    (L-l-q)  f~G   (0,qa,0)  -  G   (2a,qa,0)l 
q=l 

2       R 
+  2G  (L+l)   E  (R+l-q)  Pg   (0,0,pa)  -  G   (2a,0 ,pa)] 


+  4G   E   £   (R+l-p)  (L+l-q)  [~G   (0,qa,pa)  -  G   (2a,qa,pa)~]   . 
p=l  q=l 

(256) 


121 


(b)      Energy  of   the   array   of   F  forces; 


2F2(L+1)     .    ow2 _       ,„    _      .         „_2      L 

w      =    

F  C 


+   2F    (L+l)    G33(0,0,a)    +    2F        E    (L+l-q)[2G      (0,qa,0) 

3  "*"  q=l 


+    2G33(0,qa,a)    -   G33(0,qa,Ra)    -   G      (0  ,qa,Ra+2a) 

-   2G33(0,qa,Ra+a)J    -    (L+l)  QsG^O , 0 ,Ra+a)    +    G33(0,0,Ra) 


f    G33(0,0,Ra+2a)"l   + 


2F   ^R+1)    +    2F2(R+1)    G22(0,a,0) 


+   2F2      E   (R=l-P)  Q2G22(0,0,pa)    +    2G22(0,a,pa) 
P=l 

-G22(0,La,pa)    -  G22(0,La+2a,pa)    -   2G22(0,La+a,parj 

-  (R+l)  r"2G22(0,La+a,0)    +    G22(0,La,0)    +   G22(0 ,La+2a,ofj 

L        R  R 

-  16      E        E     G_„(0,qa,pa)    -   8      E     G      (0,La+a,pa) 

q=l   p=l  p=l 


-   8     E     G      (0,qa,Ra+a)    -   4G      (0 ,La+a, Ra+a) ; 
q=l 


(257) 


(c)   Interaction  energy  between  both  arrays: 


FG 


E   (L+l)  G   (a,0,pa)  +   E  (R+l)  Gl9(a,qa,0) 

-<        I"'  1       1^ 

Lp=l  q=l 


C  L  R 

E  (L+l-q)  Gir!(a,qa,Ra+a)  +   E  (R+l-p)  G   (a,La+a,pa) 
q=l  p=l 

-  4  ["(L+l)  G   (a, 0, Ra+a)  +  (R+a)  G   (a,La+a,0r) 


-  16 


R   L 

E   E  (L+l-q)  G   (a,qa,pa) 
Lp=l  q=l         13 


R    L 
+   E    E  (R+l-p)  G   (a,qa,pa) 
P=l  q=l 


(258) 
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The  total  energy  of  the  system  of  forces  is  the  sura  of  these  three 
equations.   After  computation  of  the  various  sums  and  collecting  terms 
containing  R  and  L,  the  following  expression  is  found  for  the  energy 
of  the  whole  array  of  forces. 


W  _  (R+l) (L+l)   (1-v)     .   ua    3.3524  -  8.6586V 


1     1    ,r      -,N  /)        2RL 


Z2  -2 


1     1     ,r,   ,%   />         2RL 


L+  >/R  +L 


/~2   2~ 
R+  VR  +L 


1     V 


"  J^~2 


2rr       2 
(1-2V) 


(L+1)  R_±VRJ±_+  (R+1)  ^LWR+L 


/2   2^ 

+  VR 


+  (R+l+L+l) 


1  (1-V)    T8.0325  -  5.5644V 


4rr 


(1-2V) 


r«.03Z5    -    5.  5644V -i 
2L  1-V  J 


1         v2 
2TT   (1-2V)2 

4TT(ia        7.  4947    -   9.  1180V 

L  s          x  - v     J 

Z2    2 

a/R   +L               1 

2                               >— 

5V--5V+1        R+L     /2     2 

1  V 


(1-2V)' 


(2.2148) 
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-(1-V)  2"      (1-V)(1-2V)2     RL 


_L   Cl-V)         f.2988  +    5.0463V-1       1  V  r5.5876    -   7.1780V~i 

+   ~^*L  ^  J+"(1-2V)2L  1'"  J 


3  V 

77    (1-2V) 


2    ' 


(259) 


If  the  array  is  extended  in  the  x  direction  such  that  L  is 
much  larger  than  R,  the  system  becomes  two  antiparallel  straight  edge 
dislocations  separated  by  a  distance  Ra.   The  energy  per  unit  length 
of  this  system  is  the  limit  of  W/(L+l)a  when  L  becomes  infinite: 
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o  2 

ua       (1-V) 


(I^l)a  -     4^     (1_2v)2 


^«  + 


(R+D 


4rrua        3.3524   -   8.6586V 
~C~  +  1-V 


H,a2    (1-V)         8.0325    -   5.5644V 


4tt 


(1-2V) 


lia  V  r  4rr(ia   _  7.4947    -  9.  1180v~j       ua_  V 


2      , 
|ia        1 

TT       1-V 


(2.2148) 


(260) 


In  the  special   case  of   V  =   1/3,    we  have 

2  *»r£c 

i 

2 


W 
(1*1)  a 


„H.a«)[^+.6»8|     -Jfc.&fc. 


2TT 


4rrixa 


-    14.9274 


(261) 


Both  force  constants  C  and  C   are  unknown  unless  C   is  taken  to  be  the 
1       o  J- 

value  found  from  the  displacement  field  computations.  However,  C3  will 
still  remain  unknown,  and  the  same  difficulty  as  was  encountered  in  the 
simulation  from  shear  loops  is  encountered  here. 

Furthermore,  Equation  (260)  does  not  have  the  well-known  form  of 
the  energy  of  such  a  system.   This  is  due  to  the  excess  strain  energy 
stored  between  the  planes  x  =  ±  a/2  after  deformation.   All  the  mate- 
rial contained  between  the  planes  x  =  ±  a  before  the  introduction  of 
the  defect  is  compressed  into  a  slab  less  twice  the  original  width. 
Unfortunately,  it  has  not  been  possible  to  propose  a  displacement  field 
in  this  region  which  would  permit  one  to  obtain  the  correction  energy 
for  this  system.   The  problem  is  much  less  simple  than  for  the  simula- 
tion from  shear  loops.   So  further  research  has  to  be  done  in  this  direction. 


CHAPTER  7 


CONCLUSIONS 


The  new  approach  presented  here  has  given  us  very  satisfying 
results  for  the  computations  of  the  small  range  displacement  field  and 
the  self-energy  of  the  various  dislocations,  especially  for  the  case  of 
the  screw  dislocation  where  a  complete  treatment  could  be  done.   For  the 
long  range  displacement  field,  this  model  is  in  agreement  with  ordinary 
continuum  model  in  all  cases.   When  compared  to  Peierls'  model,  the 
screw  dislocation  has  a  narrower  width  in  our  model,  although  the  edge 
dislocation  has  a  broader  width. 

The  arrangement  of  atoms  in  the  vicinity  of  the  dislocation  line 
is  obtained  for  most  of  the  cases  and  at  distances  closer  to  the  disloca- 
tion line  than  the  ordinary  continuum  model  permits.   It  has  not  been 
possible  to  compare  the  atomic  displacements  with  atomistic  computations, 
since  simple  cubic  structure  is  hypothetical.   Its  study,  however,  will 
be  the  guide  of  studies  on  real  crystals. 

At  this  time,  there  are  two  principal  difficulties  in  applying 
this  model:   (a)  the  evaluation  of  the  force  constant  corresponding  to 
the  forces  normal  to  the  slip  plane,  and  (b)  the  computation  of  the 
excess  strain  energy  between  the  planes  of  forces.   One  way  to  solve 
these  problems  would  be  to  compare  with  computer  simulation  models. 
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Unfortunately,  this  can  only  be  done  on  real  crystal  structures.   So  the 
next  step  for  pursuing  this  study  would  be  a  direct  application  of  the 
method  employed  here  to  reach  the  properties  of  dislocation  and  kinks 
in  face  centered  cubic  crystal,  body  centered  cubic  crystals  or  hexag- 
onal compact  crystals. 
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